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In  the  past  nine  months  we  have  continued  our  work  on  polarization  difference  vision  and 
imaging  and  their  engineering  potential  applications  to  underwater  observation  and  imaging  of 
targets  in  scattering  media.  We  have  investigated  various  aspects  of  polarization-difference 
imaging  (PDI)  system  and  its  features  and  characteristics.  Moreover,  we  have  studied 
polarization-sensitive  vision  in  certain  vertebrates,  namely,  green  sunfish  ( Lepomis  cyanellus) 
and  pigeons  ( Colwnba  livia).  In  the  following  sections,  we  report  our  Endings  in  the  past  nine 
months  in  three  projects  related  to  polarization  sensitive  vision  and  imaging. 


Reverse  Engineering  of  Polarization-Contrast  Vision: 
Polarization-Difference  Imaging  in  Scattering  Media 

In  the  past  several  months,  the  polarization  difference  imaging  (PDI)  project  has  had 
great  success  and  taken  off  in  exciting  directions.  One  phase  of  the  project  was  completed,  a 
manuscript  was  prepared,  and  new  experiments  were  designed  and  executed,  all  of  which 
showed  that  PDI  is  a  potentially  powerful  tool  for  observation  of  targets  in  scattering  media. 
This  part  of  the  report  will  cover  four  main  areas,  viz.,  observation  of  metallic  targets  in 
scattering  media,  the  development  of  an  automatic  rotational  polarizer  which  improved  the 
performance  of  the  PDI  system,  the  new  experiments  which  are  currently  underway,  and  the 
theoretical  work  which  is  proposed  for  the  near  future. 

UL _ Metallic  Targets 

At  the  end  of  last  summer,  we  were  testing  the  performance  of  the  PDI  system  on  several 
different  types  of  targets.  Initially,  we  had  the  idea  of  studying  dielectric  materials,  such  as 
plastics.  When  these  targets  are  oriented  at  an  angle  with  respect  to  the  viewing  direction  they 
give  rise  to  a  preferential  direction  of  polarization  due  to  Brewster’s  angle  phenomenon.  This 
preferential  polarization  is  usually  difficult  to  pick  up  using  conventional  polarization-blind 
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systems,  but  we  were  able  to  use  the  PDI  system  to  help  detect  these  targets  in  scattering  media 
when  the  performance  of  standard  imaging  techniques  is  pom1.  We  had  very  good  results  when 
we  tried  to  detect  transparent  targets,  and  we  had  limited  success  improving  the  detectability  of 
translucent  plastics. 

After  our  initial  experiments  on  dielectrics,  we  tried  using  the  PDI  system  on  metallic 
targets.  When  we  used  the  PDI  system  to  view  dielectrics,  we  saw  the  light  which  passed  out 
from  the  interior  of  the  target.  This  light  is  called  the  diffuse  component.  Metallic  targets  are 
different,  in  that  we  see  primarily  the  reflected  light  which  virtually  never  makes  it  into  the 
material.  Each  of  these  mechanisms  (reflection  and  diffusion)  serves  to  create  a  partial 
polarization  in  a  different  way,  but  once  the  light  leaves  the  target,  it  is  no  longer  crucial  how  the 
light  was  polarized,  only  that  it  is  and  whether  it  can  be  detected  by  our  detector.  We  first  found 
that  metallic  targets  also  partially  polarize  the  light,  but  they  were  so  visible  under  the  conditions 
of  our  experiments  that  the  PDI  system  was  not  providing  an  additional  advantage  to  improve  the 
detectability  of  the  target  relative  to  the  background.  However  upon  further  experimentations  we 
found  a  very  interesting  property  of  the  PDI  system  in  detecting  surface  features  of  metallic 
targets: 


As  an  experiment,  we  tried  scratching  the  surface  of  the  metallic  target  to  see  what  the 
effects  were.  What  started  out  as  something  of  an  unknown  turned  into  a  huge  success.  The 
scratched  surfaces  of  the  targets  acted  as  small  facets,  each  of  which  was  oriented  in  a  slightly 
different  way.  All  of  these  facets  together  gave  rise  to  a  collective  polarization  which  we  were 
able  to  detect  using  the  PDI  system.  Interestingly,  these  scratches  were  practically  impossible  to 
detect  using  standard  polarization-blind  imaging  systems  and  procedures.  However,  using  the 
PDI  system  we  were  able  to  detect  them.  Appendix  A  contains  a  preprint  of  our  manuscript 
describing  our  imaging  system  and  experimental  setup  along  with  results  showing  the  advantage 
of  polarization-difference  imaging.  Figure  1  presents  sketches  of  experimental  setup  for  such 
targets.  The  striking  outcome  of  the  PDI  system  observing  such  metallic  targets  is  shown  in  Fig. 
2  of  the  preprint  attached  in  Appendix  A. 


Fig.  1.  Top  panel  depicts  the  top  view  of  the  experimental  setup  for  our  PDI  System.  A: 
Polarization  analyzer;  TNLC:  twisted  nematic  liquid  crystal;  F:  narrow-band  filter. 
Bottom  panel  shows  the  front  view  of  the  tank  (drawn  to  scale)  with  inside  dimensions 
30cm  x  30cm  x  15cm.  This  tank  is  filled  with  water  to  which  milk  is  added.  The  target 
is  suspended  in  the  middle  of  the  tank.  The  experimental  results  of  PDI  system  in  shown 
in  Fig.  2  of  the  preprint  in  Appendix  A. 
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1.2.  Automatic  Rotational  Polarizer  (ARP) 

As  his  Senior  Undergraduate  EE  Design  Project,  Mr.  J.  Scott  Tyo  (one  of  the  student 
members  of  our  research  team)  designed  and  built  an  Automatic  Rotational  Polarizer  which 
improved  the  performance  of  our  PDI  system.  We  noticed  two  problems  which  came  up 
consistently  in  the  PDI  system  which  we  decided  to  address.  The  first  problem  was  a  motion 
artifact  that  was  due  to  the  rotating  of  our  linear  polarizer.  This  artifact  was  manifested  as  dark 
and  light  bands  around  the  edges  of  the  metallic  target.  The  second  problem  was  that  the 
background  light  may  not  be  randomly  polarized  as  it  was  in  our  early  experiments.  If  the 
background  polarization  is  not  compensated  for,  the  performance  of  the  PDI  system  can  be 
diminished. 

The  first  problem  was  solved  using  a  twisted  nematic  liquid  crystal  (TNLC)  cell.  This 
cell  serves  as  a  polarization  rotator  and  can  be  controlled  electronically.  When  there  is  a  voltage 
applied  across  the  cell,  it  does  not  affect  the  polarization  of  the  incident  light,  however,  when  the 
cell  is  not  electrically  excited,  the  polarization  of  the  light  leaving  the  cell  is  rotated  ninety 
degrees  with  respect  to  the  incident  light  We  can  take  advantage  of  the  two  states  of  the  TNLC 
to  capture  out  vertically  and  horizontally  polarized  images  without  mechanically  manipulating 
the  system.  The  introduction  of  the  TNLC  solved  our  motion  artifact  problem  and  it  was  used  to 
prepare  the  images  which  were  presented  in  Fig.  2  of  the  preprint  attached  in  Appendix  A. 

In  order  to  compensate  for  the  background  polarization  problem,  a  stepper  motor  was 
introduced.  The  TNLC  allowed  us  to  capture  images  polarized  parallel  and  perpendicular  to  our 
analyzer,  and  the  stepper  motor  rotated  our  analyzer,  thus  rotating  the  direction  of  polarization 
being  used  to  create  PD  images.  When  the  background  light  is  not  randomly  polarized,  its 
preferential  polarization  can  be  represented  as  a  vector  pointing  in  a  given  direction  (This  is  in 
fact  the  direction  of  the  major  axis  of  the  polarization  ellipse  of  polarized  part  of  the  background 
light).  In  order  to  null  out  this  background  polarization  effect  in  the  PDI  system,  the  analyzer 
must  be  oriented  so  that  it  makes  a  45  degree  angle  with  respect  to  the  major  axis  of  polarization 
ellipse  of  the  background.  The  stepper  motor  assembly  and  circuitry  has  been  designed,  built, 
and  tested,  but  it  has  not  been  used  yet  in  creating  images  for  presentation. 

L2, _ Limit  of  Detectability;  Application  of  Signal  Detection  Theory 

The  PDI  system  has  been  shown  to  improve  the  detectability  of  target  features  in 
scattering  media  (See  Fig.  2  of  Appendix  A).  We  have  observed  such  improvement 
qualitatively,  and  we  need  to  develop  some  quantitative  way  of  describing  the  advantages  of  PDI 
in  enhancement  of  target  detection  under  certain  circumstances.  In  the  areas  of  signal 
processing,  communications,  and  psychophysics,  there  is  often  the  need  to  quantify  the 
performance  of  a  given  receiver.  There  has  been  a  lot  of  work  reported  in  the  literature  in  the 
area  of  signal  detection  theory  (SDT)  [For  example,  David  M.  Green  and  John  A.  Swets,  Signal 
Detection  Theory  and  Psychophysics ,  John  Wiley  and  Sons,  Inc.,  NY  1966.].  SDT  breaks  a 
given  target  (or  stimulus  or  signal)  detection  process  down  into  the  probabilities  for  various 
events  to  occur,  such  as  hits,  false  alarms,  misses,  and  correct  rejections.  Obviously,  the 
performance  of  a  given  system  can  be  quantified,  and  in  this  way  different  systems  can  be 
compared.  In  addition,  an  ideal  observer  model  is  usually  employed  [For  example  the  work  by 
W.  S.  Geisler,  "Sequential  Ideal-Observer  Analysis  of  Visual  Discriminations"  Psychological 
Review ,  vol.  96  No.  2, 1989,  pp.  267-314.]  which  takes  into  account  only  the  information  which 
is  received,  not  any  post-receiver  operations  which  may  occur  to  degrade  a  given  signal’s 
detectability. 

We  are  in  the  process  of  applying  the  SDT  approach  to  test  the  following  hypothesis 
about  our  PDI  system:  PDI  system  can  have  a  higher  probability  of  detecting  a  target  in  a 
scattering  medium  than  standard  imaging  systems  under  certain  conditions.  These  conditions 
include  many  different  situations  where  the  signal  is  just  at  or  below  threshold  levels. 
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Obviously,  if  standard  imaging  system  can  easily  detect  a  signal,  the  PDI  system  is  not  of  high 
priority.  However  we  would  like  to  show  that  that  are  cases  where  standard  imaging  procedures 
do  not  perform  as  well  as  PDI  system  does,  and  these  cases  may  be  the  ones  with  situations  at  the 
threshold  of  detectability.  These  conditions  may  include,  but  are  not  limited  to,  environments 
where  scattering  particles  are  present  in  the  viewing  medium,  and  limit  an  imaging  system’s 
performance  and  lighting  patterns  that  have  significant  background  variations.  These  conditions 
occur  in  a  variety  of  environments,  but  two  of  the  most  common  are  underwater  and  foggy 
environments.  It  must  be  noted  that  such  possible  enhanced  detectability  for  PDI  system  appears 
to  be  in  addition  to  the  common-mode  rejection  capability  of  PD  imaging.  The  common-mode 
rejection  can  even  further  improve  the  detectability  of  targets  in  scattering  media.  (See  Fig.  2  in 
Appendix  A) 

In  order  to  test  this  hypothesis,  the  following  experiment  has  been  proposed  and  is  in  the 
early  stages  of  implementation.  Using  the  setup  described  in  Appendix  A,  we  choose  a  given 
milk  concentration.  Without  changing  anything  about  the  system,  we  collect  a  series  of 
vertically  and  horizontally  polarized  images.  From  these  images,  create  the  appropriate 
polarization-sum  (PS)  and  polarization-difference  (PD)  images.  We  will  repeat  the  procedure 
with  a  second  aluminum  target  which  is  sandblasted  but  does  not  have  the  scratched  patches. 
The  first  set  of  data  represents  the  signal  plus  noise  distribution  while  the  second  represents  the 
noise  alone  (or  non-target)  distribution.  Once  all  the  data  is  collected,  we  will  study  the 
distributions  of  various  variables,  such  as  absolute  intensity,  difference  between  different 
regions,  difference  from  the  mean,  etc.,  for  each  case.  By  plotting  these  distributions,  SDT  can 
be  used  to  analyze  the  data  and  test  the  hypothesis. 

1A _ Proposed  Theoretical  Work  for  PDI  System 

There  is  a  large  body  of  literature  which  deals  with  propagation  and  radiation  in 
scattering  and  random  media  [For  example,  A.  Ishimaru,  Wave  Propagation  and  Scattering  in 
Random  Media,  vols.  1  and  2,  Academic  Press,  NY,  1978;  Y.  Kuga  and  A.  Ishimaru,  "Imaging 
of  an  Object  Behind  Randomly  Distributed  Particles  Using  Coherent  Illumination,"  J.  Wave 
Material  Interaction,  Vol.  3,  No.  2,  pp.  105-112,  April  1988;  and  H.-W.  Chang  and  A.  Ishimaru, 
"Beam  Wave  Propagation  and  Scattering  in  Random  Media  Based  on  the  Radiative  Transfer 
Theory,"  J.  Wave-Material  Interaction,  Vol.  2,  No.  1,  pp.  41-69,  January  1987].  Some  of  these 
methods  may  be  employed  to  study  the  physics  of  the  PDI-system  experiments  presented  here. 
Specifically,  we  must  study  the  depolarization  of  a  polarized  ray  of  light  and  the  selective 
polarization  of  a  randomly  polarized  ray  of  light  Both  of  these  phenomena  are  due  to,  amongst 
other  things,  directional  scattering  of  light  within  the  scattering  medium.  We  would  like  to  study 
ways  of  analytically  justifying  the  results  which  we  obtain  from  the  SDT  experiments.  This 
justification  will  probably  come  out  of  the  scattering  medium  work,  but  new  methods  may  also 
be  employed. 


Biophysical  Basis  of  Polarization  Sensitivity  in  Teleost  Fish 

The  goal  of  this  part  of  the  research  program  is  to  test  the  hypothesis  that  sunfish  double 
cone  photoreceptors  behave  as  birefringent  waveguides,  and  thus  provide  the  mechanism  by 
which  these  fish  detect  polarization  of  light.  In  previous  work,  (See  Appendix  B  for  a  reprint  of 
our  previous  work)  electron  microscopic  analysis  of  sunfish  cone  inner  segments  led  to  the 
hypothesis  that  the  inner  segment  region  of  the  receptors  are  graded  index  waveguides.  During 
this  reporting  period  we  collected  data  to  test  that  hypothesis  as  directly  as  possible  given 
currently  available  technology. 

One  of  us  (Mickey  Rowe  who  is  one  of  our  graduate  students  from  The  Institute  of 
Neurological  Science  at  the  University  of  Pennsylvania)  made  three  trips  to  the  Duke  University 
Medical  Center  in  Durham,  North  Carolina  to  collect  interferometric  data  from  sunfish  cones. 


( 
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He  collected  the  data  in  the  laboratory  of  Dr.  Joseph  Corless,  whose  modified  Vickers  M8S 
microscope  made  the  data  collection  possible. 

2.1  Scanning  Interferometry 

After  the  sunfish's  retina  is  obtained,  it  is  left  in  a  Ringer’s  solution  similar  to  that  used  by 
Dearry  and  Burnside  [A.  Deany  and  B.  Burnside,  "Effects  of  Extracellular  Ca++,  K+,  and  Na+ 
on  Cone  Retinal  Pigment  Epithelium  Retinomotor  Movements  in  Isolated  Teleost  Retinas",  J. 
General  Physiology,  Vol.  83,  pp.  589-611, 1984].  It  was  discovered  during  the  third  trip  that  the 
best  results  were  obtained  by  shaking  the  retina  in  the  solution,  and  pipetting  a  small  volume  (8 

pi)  of  fluid  from  the  surrounding  solution.  It  was  also  discovered  just  prior  to  the  third  trip  that 
the  largest  number  of  cones  could  be  obtained  by  killing  the  fish  2-3  hours  into  its  subjective 
night  alter  the  fish  had  been  in  complete  darkness  for  at  least  26  hours.  During  the  first  two  trips 
very  little  data  were  collected  due  mainly  to  problems  with  the  dissection. 

The  Vickers  M85  microscope's  operation  can  be  briefly  described  as  follows:  light  from 
a  He-Ne  laser  (wavelength  632.8  nm  in  free  space)  is  split  into  two  beams,  a  scanning  beam  and 
a  reference  beam.  The  beams  travel  in  parallel  through  the  microscope’s  optical  system  until 
after  they  have  passed  through  the  sample  plane  in  which  the  photoreceptors  reside.  While 
travelling  through  the  sample  plane,  one  of  the  beams  passes  through  a  photoreceptor  of  interest, 
and  the  other  beam  passes  through  the  solution  adjacent  to  the  receptor.  Passage  through  the 
photoreceptor  retards  the  scanning  beam,  so  that  when  the  two  beams  are  recombined  beneath 
the  sample  plane,  the  scanning  beam  is  shifted  in  phase  relative  to  the  reference  beam.  The 
microscope  is  designed  to  measure  this  phase  shift  and  thus  determine  the  optical  path  length  of 
the  section  of  the  photoreceptor  traversed  by  the  beam.  For  these  experiments,  the  beam  was 

reduced  to  a  spot  size  of  approximately  0.3  pm  as  it  passed  through  the  photoreceptor.  During  an 

experiment  the  beams  are  moved  in  a  raster  pattern  in  the  sample  plane,  with  0.22  pm  between 
measurements.  Data  are  thus  collected  as  a  two  dimensional  matrix  of  points  (see  Fig.  2). 

The  data  collected  as  described  in  the  preceding  paragraph  represent  the  additional  optical 
path  length  (opl)  travelled  by  the  scanning  beam  relative  to  that  travelled  by  the  reference  beam. 

This  opl  is  related  to  the  refractive  index  via:  opl(x,z)  =  J 5 n(x,y,z)dy.1  (See  Fig.  2  for  an 

j 

explanation  of  the  coordinate  system  used  throughout  this  section.)  For  this  project  we  are 
interested  in  extracting  n(x,y,z)  from  our  data. 


Although  our  goal  is  to  understand  the  biophysics  of  the  double  cones,  we  have  focused 
our  efforts  on  analyzing  data  collected  from  single  cones  because  they  can  be  assumed  to  be 
radially  symmetric.  This  assumption  of  radial  symmetry  allows  us  potentially  to  extract  n(x,y,z) 
from  a  single  matrix  of  opl(x,z)  data.  The  previous  evidence  which  lead  us  to  propose  a  gradient 
index  model  of  the  double  cones  applies  equally  well  to  the  single  cones.  If  it  can  be  shown  that 
the  single  cones  possess  an  index  profile  graded  as  we  have  hypothesized,  then  we  can  infer  that 
the  double  cones  have  a  graded  index  profile  as  well. 


1  The  6  in  this  and  subsequent  equations  signifies  a  difference  between  the  refractive  index  of  the  object  under 
study,  and  the  Ringer’s  solution  in  which  the  object  is  immersed.  The  Ringer’s  solution  had  a  measured  refractive 
index  of  1.3351,  aid  this  constant  can  be  added  to  any  of  the  Sn's  to  convert  them  to  absolute  refractive  indices. 
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Figure  2.  Diagramatic  representation  of  scanning  interferometry.  The  scanning  and 
reference  beams  arc  moved  in  a  raster  pattern  in  order  to  collect  a  two-dimensional  matrix 
of  data  samples.  At  each  data  point,  the  phase  difference  between  the  scanning  and 
reference  beams  indicates  the  additional  optical  path  traveled  by  the  scanning  beam  as  it 
travels  through  die  cone  or  other  specimen  of  interest 


2.2  Prelir, 


y  Results 


In  our  data  analysis  we  have  made  some  simplifying  assumptions,  chiefly  that  the 
scanning  beam  passes  through  our  samples  as  diagrammed  in  Fig.  2.  In  reality,  the  scanning 
beam  is  too  narrow  to  be  treated  as  a  pencil  of  light— because  the  spot  is  focused  to  a  point 
smaller  than  the  light’s  wavelength,  the  “beam”  must  be  diverging  at  the  sample  plane.  The 
microscope  is  designed  well  enough  (e.g.  with  appropriately  placed  apertures)  that  this  may  not 
affect  the  calculations  presented  below,  but  we  are  in  the  process  of  performing  some  additional 
calibrations  to  test  the  validity  of  this  assumption.  Even  if  the  light  hitting  the  sample  was  a 
straight  beam,  refraction  at  the  interface  should  cause  the  beam  to  bend  as  it  enters  the  cones  at 
all  positions  except  the  midline  (x  =  0,  where  the  light  impinges  normally  on  the  surface).  With 
preliminary  analyses  based  on  a  homogeneous  refractive  index  model  this  bending  was  found  to 
be  minor,  and  it  is  thus  reasonable  to  ignore  it  at  this  stage  of  calculation.  Furthermore,  the 
increased  pathlength  at  the  periphery  of  the  cell  would  cause  that  region  to  appear  to  have  a 
larger  refractive  index  than  it  does,  and  thus  would  bias  our  results  away  from  the  conclusion  we 
are  attempting  to  support. 
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Given  the  assumptions  that  the  single  cones  are  radially  symmetric  around  their  z-axes, 
and  that  the  scanning  beam  passes  through  them  as  an  undeflected  pencil  of  light,  we  can  extract 
n(x,y,z)  from  opl(x,z)  in  at  least  two  different  ways.  The  way  that  we  have  concentrated  upon 

during  this  reporting  period  is  to  begin  with  models  of  n(r,z)  (where  r  =  ^x2  +  y2 ),  determine 
the  theoretical  opl(x,z)  which  would  result,  and  calculate  the  difference  between  the  data  and  this 
model  ami  square  that  difference.  We  then  adjusted  parameters  of  the  model  and  found  the 
squared  error  again.  This  procedure  was  repeated  until  a  minimum  in  the  difference  between 
data  and  theory  were  obtained.  Prior  to  fitting,  however,  the  data  were  made  symmetrical.  Since 
we  are  assuming  radial  symmetry,  it  should  be  true  that:  opl(x,z)  =  opl(-x,z).  We  thus  treated 
the  data  on  either  side  of  the  z  axis  as  independent  estimates  of  the  actual  value.  In  order  to 
appropriately  average  the  data,  we  had  to  determine  where  x  =  0.  We  did  this  independently  for 

each  value  of  z  in  our  data  set  by  maximizing  j"  opl(x)  •  opl(-x)dx  as  the  position  where 

x  =  0  was  moved  across  the  data  slice.  Since  our  data  is  discrete  and  the  position  where  x  =  0 
might  not  fall  on  a  pixel,  data  were  interpolated  prior  to  averaging.  The  resulting  “symmetrized” 
data  were  used  in  all  subsequent  calculations.  The  first  index  model  tested  was: 


homogeneous:  8n(x:  z)  = 

0  else 

The  homogeneous  model  has  two  parameters,  Snmax  and  R»  the  refractive  index  within,  and  the 
radius  of  the  cell  respectively.  Note  that  both  parameters  are  potentially  functions  of  z.  The 
second  model  tested  was: 

parabolic:  5n(x:z)  =  “  (^r)  ]  lxl  ^  R 

0  else 


This  model  also  has  two  parameters,  but  in  this  case  Snmax  is  the  refractive  index  increment  at 
the  center  of  the  cell  and  R  specifies  not  only  the  outer  border  of  the  cell,  but  also  how  fast  the 
index  declines  as  a  function  of  r. 


Overall,  the  data  were  better  fit  by  the  homogeneous  model  than  the  parabolic  model. 
However,  when  the  models  were  compared  point  by  point  it  was  determined  that  around  x  =  0, 
the  parabolic  model  over-estimated  the  opl  and  the  homogeneous  model  under-estimated  it. 
Nearer  to  the  perimeter  of  the  cells  the  situation  was  reversed:  the  parabolic  model  tended  to 
under-estimate  the  opl  while  the  homogeneous  model  over-estimated  it.  For  this  reason  we 
chose  to  pursue  a  three-parameter  model: 


power-law: 


Sn(x:z)  =  6n“>[l-(x/R)”]  MSR 

0  else 


The  power-law  model  is  a  more  general  expression  of  the  two  particular  models  already 
described.  The  power-law  model  reduces  to  the  homogeneous  model  when  a  =  <»,  and  to  the 

parabolic  model  when  a  =  2.  For  intermediate  values  of  a,  the  index  profile  has  a  shape 
intermediate  between  the  two  previous  cases.  After  fitting  the  power-law  model  to  the  data,  the 
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values  for  naVg,2,R,  and  a  were  plotted  as  a  function  of  z.  Figure  3  ^nows  the  results  of  these 
calculations  along  with  the  results  from  the  homogeneous  model  for  one  cell. 

The  programs  currently  used  to  find  the  best  fits  for  the  power-law  model  run  very 
slowly.  The  cell  depicted  is  currently  the  only  one  for  which  the  analysis  has  been  completed. 
However,  for  the  homogeneous  model,  several  other  cells  have  been  tested.  The  general  trend 
displayed  in  panel  A  of  Fig.  3  is  evident  in  nearly  all  of  these  cells.  That  is,  the  refractive  index 
is  tower  near  the  base  of  the  photoreceptor  inner  segment  than  it  is  near  die  base  of  the  outer 
segment,  and  it  rises  monotonically  in  between.  This  pattern  provides  continuation  of  a 
prediction  made  previously  based  on  electron  microscopy  alone  (Our  paper  in  Appendix  B).  The 
sharp  drop  in  index  at  the  outer  segment  is  unexpected,  and  may  be  the  result  of  post-mortem 
pathology;  the  outer  segment  is  more  fragile  and  thus  more  sensitive  to  trauma  imposed  by  the 
dissection  procedure  than  is  the  inner  segment.  Panel  B  depicts  the  dependence  of  the  cell’s 
radius  on  z.  The  values  obtained  via  this  method  are  similar  to  light  microscopical 
measurements  of  fixed  and  stained  retinas  (D.  A.  Cameron  and  E.  N.  Pugh,  Jr.,  Nature ,  vol.  353, 
161-164, 1991;  and  D.  A.  Cameron  and  S.  S.  Easter,  Jr.,  Vision  Neuroscience,  vol.  10,  375-384, 
1993).  Panel  C  of  Fig.  3  provides  support  for  the  prediction  we  set  out  to  test.  Two  anticipated 
points  are  demonstrated  in  that  panel:  1)  the  cell  is  best  tit  by  values  of  alpha  less  than  infinity. 
If  the  cell  really  had  a  homogeneous  refractive  index  then  one  would  predict  that  no  minimum 
error  would  be  found— the  error  in  the  best  tit  would  continue  to  decline  as  higher  values  of  alpha 
were  tested.  2)  The  best  fit  alpha  increases  as  a  function  of  z.  This  is  in  accord  with  previous 
observations  based  on  electron  microscopy. 

For  a  quantitative  evaluation  of  the  tits  for  the  homogeneous  versus  power  law  models, 
the  error  functions  were  plotted  as  a  function  of  z.  As  described  above,  the  error  between  the 
data  and  the  model  were  determined  by  finding  the  opl  which  would  result  from  a  model  having 
the  parameters  depicted  in  Fig.  3,  subtracting  those  values  from  the  symmetrized  opl  data,  and 
then  squaring  the  resultant  number  at  each  point  (x,z).  These  numbers  were  then  summed  along 
x  to  show  the  error  only  as  a  function  of  z.  Fig.  4  portrays  the  residual  errors  after  the  best  fits 
were  obtained  for  the  two  types  of  model.  Note  that  the  errors  arc  plotted  on  a  logarithmic  scale. 
The  errors  have  not  been  normalized-the  errors  seem  to  climb  as  a  function  of  z,  but  this  may 
result  from  the  fact  that  the  raw  numbers  in  the  data  increase  as  a  function  of  z.  However,  the 
important  point  to  note  is  that  the  power-law  model  tends  to  fit  the  data  around  ten  times  as  well 
as  the  homogeneous  model  along  most  of  the  photoreceptor  inner  segment's  length. 

Another  method  for  extracting  n(r,z)  from  opl(x,z)  is  to  perform  an  inverse  Abel 
transform  of  the  data.  The  Abel  transform  is  a  special  case  of  the  Radon  transform  used  in 

medical  imaging.  It  can  be  shown  that  if  8n  is  radially  symmetric,  it  can  be  reconstructed 
according  to  Barrett  and  Swindell3 : 


8n(r)  = 


where  in  our  calculation  the  upper  limit  is  approximately  taken  to  be  R  which  should  be  obtained 
from  the  data.  The  difficulty  in  applying  this  method  is  that  d(opl)/dx,  and  R  must  be  extracted 
from  the  data.  For  one  scan  line  of  the  same  cell  whose  analysis  is  depicted  in  Figs.  3-4,  we  have 


2  Snavg  is  related  to  8nmax  in  the  power-law  model  via:  8nml  =  +  2)  *  This  08,1  **  shown  by 

integrating  8n(x)  from  -R  toR,  and  dividing  by  xR2. 

3  H.  H.  Barrett  and  W.  Swindell,  Radiological  Imaging;  The  Theory  of  Image  Formation,  Detection,  and 
Processing,  vol.  2,  Academic  Press,  NY,  1981,  p.  406. 
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Figure  3.  Best  fit  parameters  for  the  homogeneous  and  power  law  models  of  the 
refractive  index  of  a  single  cone.  A.  Average  refractive  index  within  the  single  cone  as  a 
function  of  distance  from  the  base  of  the  inner  segment  The  large  drop  near  z  =  14pm  is 
the  junction  between  the  inner  and  outer  segment  of  the  photoreceptor.  B.  Radius  of  the 
cone.  C.  For  the  power-law  model,  the  exponent  specifying  how  fast  n  decreases. 

performed  this  calculation  by  fitting  the  opl  data  to  a  spline  curve,  calculating  the  derivative  from 
this  spline. 


gradient  model 
homogeneous  model 


Figure  4.  Residual  error  between  models  and  data.  For  each  point  along  the  z-axis,  optical 
path  lengths  were  calculated  as  a  function  of  x  to  best  fit  the  data.  The  difference  between 
the  calculated  opl’s  and  the  symmetrized  data  were  then  found  as  a  function  of  x.  The 
differences  were  then  squared  point  by  point  and  the  resultant  squared  error  summed  over  x. 
These  values  are  here  plotted  as  a  function  of  z  for  the  homogeneous  and  power-law 
models. 
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and  then  defining  R  as  the  appropriate  zero  crossing  of  the  d(opl)/dx  curve.  The  resulting  plot  of 
n(r)  is  depicted  in  Fig.  5  with  the  n(r)  determined  from  the  power  law-model  fit  included  for 
comparison. 

At  this  time  it  is  not  clear  whether  the  small  central  rise  in  n(r)  as  determined  by  the 
inverse  Abel  method  is  an  accurate  reflection  of  the  underlying  refractive  index  profile,  or 
whether  it  is  an  artifact  of  the  reconstruction  technique.  One  line  of  evidence  suggesting  that  the 
hump  is  real  (and  therefore  that  the  power-law  function  we  are  using  is  incapable  of  faithfully 
recovering  the  true  index  profile)  is  that  when  the  best  fit  opl  is  subtracted  from  the  data,  the 
model  still  seems  to  underestimate  the  center  and  overestimate  the  perimeter.  Further  analysis  of 
the  data  will  hopefully  determine  if  another  model  can  fit  the  data  better  than  the  power-law 
function. 


inverse  Abel 
gradient  fit 


x(Hm) 


Figure  5.  n(r)  was  reconstructed  from  opl(x)  for  one  scan  line  of  the  cell  analyzed  for 
Figures  2-4.  The  reconstruction  was  done  in  two  ways.  The  solid  curve  results  from 
applying  the  inverse  Abel  transform  on  the  symmetrized  data.  The  dotted  curve  is  from 
the  best  fit  with  a  power-law  model  The  scan  line  in  question  is  for  z  »  8.36  pm. 


2 J - Behavioral  Experiment  of  Polarization  Sensitivity  of  Green  Sunfish 

We  have  initiated  work  on  behavioral  experiments  studying  sensitivity  of  green  sunfish's 
visual  system  to  the  polarized  light.  An  undergraduate  student,  Ms.  Kara  Ko  has  joined  our 
group  over  the  summer  to  work,  along  with  our  graduate  students,  in  this  part  of  the  project.  We 
are  exploring  several  possibilities,  namely,  the  ability  of  fish  to  discriminate  between  linearly 
polarized  light  and  circularly  polarized  light,  to  discriminate  between  vertically  and  horizontally 
polarized  light,  and  also  to  discriminate  between  monochromatic  lights  with  differing 
wavelengths.  Experimental  setup  is  being  prepared  for  measurement  of  electroretinographic 
(ERG)  signals  from  fish's  eye  when  it  is  illuminated  by  flashes  of  light. 
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Biophysical  Basis  of  Polarization  Sensitivity  in  Pigeons 

Modifications  and  improvement  of  new  apparatus  for  pigeon  experiments  (for  lateral 
testing)  and  the  original  (frontal)  were  completed  in  September  93.  The  computer  programs 
necessary  to  run  this  new  apparatus  were  written.  Tests  on  one  bird  were  run  on  the  frontal 
contrast  paradigm.  In  this  contrast  paradigm,  two  possible  stimuli  were  the  patterned  field  (i.e., 
left  half  of  rectangular  field  was  black  and  the  right  half  was  white),  and  the  homogeneous  field 
(i.e.,  left  part  was  white,  and  the  right  part  was  also  white  with  a  black  line  separating  the  two 
fields).  The  homogenous  field  was  S+  and  the  patterned  field  was  S_.  This  paradigm  was  used 
for  two  reasons:  1)  we  needed  a  control  for  the  experiment  we  ran  last  year,  in  which  the  S+  was 
the  patterned  field  and  S_  was  homogenous  one;  and  2)  we  needed  more  subjects  run  in  this 
experiment  We  did  not  think  that  one  subject  was  enough.  The  results  were  comparable  to  that 
run  in  the  previous  year  (For  70  %  of  the  trials,  the  performance  of  the  bird  was  correct  when  the 
contrast  between  the  two  sides  of  the  test  field  was  around  21  %.) 

One  bird  was  run  on  the  lateral  testing  paradigm.  From  October  19  -November  15,  this 
bird  was  trained  to  peck  at  the  orienting  key.  Initially,  one  peck  was  resulted  in  the  food  reward. 
The  number  of  pecks  needed  was  slowly  increased.  Eventually,  we  paired  the  orienting  key  with 
S+  test  field  presentation.  So,  both  were  simultaneously  illuminated.  It  was  hoped  that  the 
animal  would  happen  to  look  down  the  side  channel  while  pecking  the  orienting  key.  Once  the 
animal  learned  this  task,  testing  was  begun.  It  turned  out  that  it  was  very  difficult  to  perform 
such  tests.  The  main  reason  was  because  we  were  not  certain  that  the  animal  was  looking  down 
the  test  channel.  In  an  attempt  to  make  the  orienting  key  serve  as  an  indicator  of  when  to  "look" 
down  the  test  channel,  we  gave  it  primary  importance  over  the  channel.  The  "strength"  of  the 
orienting  key  as  a  signal  to  when  food  was  going  to  be  presented  overrode  the  importance  of 
looking  down  the  test  channel  to  examine  the  state  of  the  test  field. 

In  February,  a  couple  of  different  methods  were  used  in  order  to  get  the  birds  to  pay  more 
attention  to  the  test  channel.  First,  we  increased  the  no-reward  option  if  the  animal  pecked 
during  the  S_  condition.  This  didn't  work.  Next,  we  decided  that  keeping  the  orienting  key 
steadily  illuminated  was  a  better  method  of  testing  the  animal.  In  this  way  the  state  of  the 
orienting  key  did  not  change,  but  the  state  of  the  channel  did.  First,  a  tone  was  added  to  signal 
the  presentation  of  the  test  field.  So,  while  the  animal  was  pecking  at  the  orienting  key,  two  250- 
ms  tone  burst  were  emitted  from  a  location  directly  above  die  test  field.  500-ms  after  the  second 
tone  burst,  the  test  field  was  turned  on.  During  the  S_  stimulus  presentation,  the  holder  which 
mounts  either  the  polarizer  or  quarter-wave  plate  is  rotated.  In  the  S+  condition,  the  field  is  held 
constant.  For  two  birds,  a  contrast  paradigm  was  used,  so  the  S+  was  the  steady  homogeneous 
(white)  field  and  the  S_  was  the  rotating  field.  For  one  bird,  the  S+  was  the  circularly  polarized 
field,  and  the  S_  was  the  rotating  polarization  field.  The  plates  rotated  at  0.6  Hz.  During  S_ 
presentation,  the  animal  was  supposed  to  inhibit  pecking.  If  the  bird  pecks  the  key  the  same 
number  of  times  as  required  in  the  S+  condition,  it  is  sent  into  a  no-reward  phase.  In  the  no¬ 
reward  phase,  a  loud  500-ms  tone  immediately  follows  each  additional  peck  in  order  to  call 
attention  to  the  lateral  field  and  illustrate  that  pecking  during  this  stimulus  state  was  unwanted. 
The  birds  began  showing  some  inhibition  in  this  paradigm,  but  the  difference  in  pecking  rates  did 
not  vary  enough  between  the  two  stimulus  conditions. 

In  the  middle  of  March  we  changed  the  design  of  the  experiment  once  again.  The  main 
difference  occurs  during  the  S_  presentation.  Basically,  we  have  added  some  training  S_  trials. 
Training  trials  are  concentrated  in  the  beginning  of  the  experimental  session,  with  a  few  scattered 
in  the  middle  and  the  end.  In  training  trials,  two  pecks  at  the  orienting  key  sends  the  animal  into 
the  no-reward  phase  described  above.  The  purpose  of  these  trials  is  to  extinguish  pecking  to  S_ 
early  in  the  experimental  session.  All  three  birds  were  tested  using  the  rotating  contrast 
paradigm  just  described  above.  The  one  bird  was  switched  over  from  the  pure  polarization 
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paradigm  because  we  thought  it  was  a  much  harder  task,  and  wanted  to  get  her  to  attend  to  the 
lateral  field.  All  three  birds  are  showing  significant  inhibition  during  S_  presentation.  The 
pecking  rates  during  S+  and  S_  conditions  are  statistically  significant. 

The  problem  with  this  experiment  is  that  the  rotating  mount  serves  as  an  audible  cue  for 
the  animal.  We  are  presently  trying  to  eliminate  this  as  a  cue.  What  we  want  to  do  is  to  have  the 
mount  rotate  in  both  conditions  because  we  think  that  a  flickering  stimulus  may  be  easier  for  the 
animal  to  detect  Recently,  we  bought  a  motor  which  enables  us  to  rotate  the  stimulus  at  a  flicker 
rate  up  to  10-Hz.  However,  we  have  been  unable  to  implement  this  new  procedure  because  of 
the  physical  features  of  the  light.  In  order  to  move  the  rotatable  mount  during  the  homogeneous 
field  presentation,  the  light  must  be  100%  circularly  polarized  before  it  reaches  the  rotating 
polarizer.  However,  due  to  optical  constraints,  we  have  not  yet  been  able  to  achieve  this.  As  a 
result,  the  stimulus  state  is  not  static  as  we  rotate  the  polarizer.  We  are  currently  attempting  to 
solve  this  problem. 
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Preprint  of  the  Manuscript  entitled 

Polarization-Difference  imaging  (PDI): 
a  biologically  Inspired  Technique  for  Observation  through 

Scattering  Media 

Authors:  M.  P.  Rowe,  J.  S.  Tyo,  E.  N.  Pugh,  Jr.,  and  N.  Engheta 
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ABSTRACT 

Many  animals  have  visual  systems  that  exploit  the  polarization  of  light,  and  some  of  these 
systems  are  thought  to  compute  difference  signals  in  parallel  from  arrays  of  photoreceptors 
optimally  tuned  to  orthogonal  polarizations.  We  have  hypothesized  that  such  "polarization- 
difference"  systems  may  improve  the  detectability  of  objects  in  scattering  media  by  reducing  the 
effectiveness  of  background  light  whose  degree  of  linear  polarization  is  different  from  that  of  the 
objects.  Here  we  present  experimental  results  which  demonstrate  that  a  man-made  polarization- 
difference  imaging  system  can  enhance  the  visibility  of  features  of  a  target  in  a  scattering 
medium.  Specifically,  we  show  that  surface  features  of  a  metallic  target  invisible  to  conventional 
imaging  are  readily  visible  with  polarization-difference  imaging. 


Optical  scattering  by  suspended  particles  (e.g.  fog,  rain,  plankton)  is  well  known  to 
diminish  the  visual  contrast  of  objects.1*3  Although  polarization-sensitive  vision  is  well 
documented  as  serving  in  navigation1*4*5,  some  types  of  polarization-sensitive  vision  may  also 
serve  to  enhance  the  visibility  of  targets  in  scattering  media.6*7  The  goal  of  this  Letter  is  to 
demonstrate  that  a  man-made  polarization-difference  imaging  (PDI)  system  similar  to  that 
hypothesized  to  function  in  some  biological  visual  systems7*8  can  indeed  enhance  the  visibility 
of  target  features  in  a  scattering  medium. 

Our  PDI  system  captures  images  of  a  given  scene  at  two  orthogonal  linear  polarizations, 
computes  the  difference  of  the  images  pixel-by-pixel,  and  re-scales  the  resultant  difference- 
image  for  optimum  use  of  the  display  range.  Symbolizing  the  two  image  intensity  distributions 
as  I||(x,y)  and  Ij_(x,y),  where  (x,y)  identifies  the  pixel  position  in  the  image  and  II  and  X  indicate 
two  orthogonal  linear  polarizations  defined  below,  the  PDI  system  generates  the  polarization- 
difference  (PD)  image 

pdl(x,y)  =  I||(x,y)-Ii(x,y).  (la) 

The  system  also  computes  the  polarization-sum  (PS)  image 

psl(x,y)  =  I|j(x,y)  +  Ix(x,y).  (lb) 

For  an  ideal  linear  polarizer,  the  PS  image  is  equivalent  to  a  polarization-blind  image  obtained  by 
a  conventional  imaging  system.  We  define  the  local  "observed  degree  of  linear  polarization"  by 

ODLP(x,y)  h  pdI(x,y)/pSI(x,y).  (2) 

The  quantity  <ODLP>region  =  <pdl(x,y)>regi0n/<psl(x  ,y)>region  serves  as  a  dimensionless 
measure  of  the  observed  PD  magnitude  in  a  particular  image  region.  In  general,  both  pdl(x,y) 
and  ODLP(x,y)  depend  upon  the  choice  of  polarization  axes;  psl(x,y),  however,  does  not. 

The  experimental  layout  is  sketched  in  Fig.  1.  Two  incandescent  tungsten  filament  slide 
projectors  back-illuminate  a  sheet  of  0.25"  thick  white  plexiglass  attached  to  one  face  of  a  glass 
tank.  The  tank  is  filled  with  water  to  which  milk  is  added.  The  plexiglass  acts  as  an  initial 
diffuser  and  the  milk  as  a  scattering  agent.  A  CCD  camera  (Hamamatsu,  model  XC-77)  views 
the  tank  through  the  clear  glass  face.  In  the  absence  of  a  target,  the  radiance  of  the  tank  is  nearly 
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uniform  across  the  viewing  area,  and  is  essentially  randomly  polarized.  The  target  used  for  the 
experiments  reported  here  is  an  aluminum  disk.  The  target  is  suspended  in  the  middle  of  the  tank 
by  means  of  a  rod  attached  to  the  disk  face  occluded  from  the  camera.  The  rod  is  attached  to  a 
clear  plexiglass  mount  suspended  from  above;  the  surface  of  this  mount  is  perpendicular  to  the 
incident  light,  and  both  it  and  the  rod  attached  to  the  disk  are  undetectable  to  both  PDI  and 
conventional  imaging.  The  disk  surface  facing  the  camera  is  sandblasted,  rendering  it  nearly 
Lambertian  everywhere,  except  for  two  1-cm2  patches.  These  patches  are  raised  a  few 
thousandths  of  an  inch  and  have  been  abraded  with  emory  paper  lighdy  in  orthogonal  directions 
--  on  one  patch  roughly  radially,  and  on  the  other  roughly  tangentially.  The  camera  is  equipped 
with  a  macro  lens  (Vivitar  Y/C  55  mm  f2.8)  and  forms  an  image  of  the  target  plane  (see  Figure 
2);  images  are  digitized  and  processed  with  an  image  analysis  system  (Series  151,  Imaging 
Technology,  Inc.).  Before  reaching  the  camera,  the  light  from  the  tank  passes  through  a  twisted 
nematic  liquid  crystal  (TNLC-Liquid  Crystal  Institute,  Kent  State  U.)  In  its  off  state  the  TNLC 
rotates  the  plane  of  polarization  of  incident  light  by  90°;  when  driven,  the  TNLC  passes  the 
incident  light  with  no  rotation.  After  the  TNLC,  the  light  passes  through  a  linear  analyzer.  Since 
the  TNLCs  behavior  is  wavelength-dependent,  a  narrow-band  filter  (10  nm  full-width  at  half 
maximum,  centered  at  610  nm)  serves  to  eliminate  light  out  of  the  TNLCs  operating  waveband. 
The  TNLC/analyzer  combination  selects  the  orthogonal  axes  of  the  polarization-difference  image 
as  follows:  when  the  TNLC  is  driven,  the  TNLC/analyzer  passes  only  light  polarized  parallel  to 
the  analyzer  axis,  yielding  the  image  IiKx.y);  when  the  TNLC  is  off,  the  TNLC/analyzer 
combination  passes  light  polarized  perpendicular  to  the  analyzer  axis,  yielding  I±(x,y).  For  the 
results  reported  here,  the  analyzer  axis  was  oriented  vertically. 

Figure  2  illustrates  the  application  of  the  PDI  system  to  the  aluminum  target  suspended  in 
diluted  milk.  Panels  A  and  B  present  the  images,  In(x,y)  and  Ii(x,y),  convolved  with  a  two- 
dimensional  (low-pass)  filter.  Panels  C  and  D  represent  the  sum  and  difference  images.  Panels 
E  and  F  present  the  data  in  images  in  C  and  D,  but  now  scaled  to  utilize  the  full  8-bit  intensity 
range  of  the  display.  The  abraded  patches,  which  are  not  visible  in  Panels  A,  B  and  C  (the  raw 
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and  sum  images)  are  clearly  visible  in  the  scaled  PD  image  (Panel  F),  but  barely  visible  in  the 
scaled  PS  image  (Panel  E).  Panels  A'  -  F  present  numerical  plots  (as  a  function  of  x)  of  average 
pixel  intensities  in  the  vertical  bands  between  the  arrows  shown  in  A  -  F;  these  plots  provide 
quantitative  evidence  supporting  the  qualitative  conclusion  drawn  from  inspection  of  the  images, 
viz.,  that  the  abraded  patch  regions  are  most  distinct  from  the  disk  background  in  the  scaled 
difference  image  (F,  F). 

The  principal  factor  underlying  the  enhanced  visibility  of  the  two  patches  in  Panel  F  is 
the  common-mode  rejection  feature  intrinsic  to  PDI.  Given  that  a  target  or  target  feature 
produces  a  non-zero  value  of  <ODLP>  in  some  region  (Eq.  2),  the  system  enables  extraction  of 
this  feature  by  "rejecting"  intensity  "common"  to  both  polarization  axes.  This  common-mode 
rejection  feature  of  the  PDI  system  is  exhibited  in  two  ways  in  the  images  of  Fig.  2:  firstly,  the 
relatively  intense  halo  of  unpolarized  light  surrounding  the  disk  is  eliminated  (cf.  C  and  D'),  and 
secondly,  the  relatively  modest  (ca.  20%)  intensity  variation  across  the  disk  surface  is  also 
minimized,  allowing  higher  gain  to  be  applied  to  the  final  displayed  PD  image  (cf.  E'  and  F). 

Based  on  the  results  shown  in  Fig.  2  and  others  we  have  obtained,  PDI  can  be  argued  to 
have  a  very  general  applicability.  Firstly,  our  results  show  that  PDI  is  quite  sensitive  to  small 
signals.  For  the  left  abraded  patch  in  Fig.  2,  <ODLP>  =  0.0164;  for  the  right  patch,  <ODLP>  = 
-0.0138.  Moreover,  in  experiments  in  which  we  systematically  raised  the  milk  concentration  to 
degrade  the  images  until  the  target  patches  were  undetectable,  a  target  patch  having  <ODLP>  < 
0.01  could  still  be  easily  detected  in  the  PD  image  (data  not  shown).  The  images  of  great  many 
object  surfaces  in  natural  environments  predictably  will  have  ODLPs  of  considerably  higher 
magnitude.9 

Secondly,  PDI  generalizes  readily  to  scattering  environments  in  which  the  "background" 
itself  has  some  degree  of  linear  polarization.  Judicious  selection  of  the  orthogonal  polarization 
axes  of  the  PDI  system  can  enhance  the  ODLP  of  target  regions  relative  to  that  of  the 
background  in  the  image  plane.  One  of  us  (J.  S.  Tyo,  EE  Senior  Design  Project)  has  produced  a 
system  that  automatically  finds  the  polarization  axes  that  minimize  the  ODLP  under  "background 
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alone"  conditions;  once  <ODLP>t,ackground  is  minimized,  maximum  display  gain  can  be 
automatically  introduced,  and  targets  having  an  ODLP  different  from  the  background  readily 
"pop  out." 

Finally,  PDI  also  possesses  the  generally  useful  qualities  of  being  passive,  simple  and 
potentially  very  fast  PDI  can  operate  passively  in  any  region  of  the  electromagnetic  spectrum 
in  which  natural  radiation  exists.  PDI  is  simple,  in  that  it  does  not  require  the  use  of 
sophisticated  image  processing  techniques  (but  of  course,  any  image  processing  techniques  can 
be  added  to,  or  used  in  parallel  with  PDI  to  make  it  even  more  effective).  PDI  is  also  potentially 
very  fast,  inasmuch  as  it  can  be  implemented  in  a  massively  parallel  system,  as  nature  appears  to 
have  done  in  the  retinas  of  many  animals.7*8*10 
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Figure  Captions 

Fig.  1.  Top  panel  depicts  the  top  view  of  the  experimental  setup.  A:  Polarization  analyzer, 
TNLC:  twisted  nematic  liquid  crystal;  F:  narrow-band  filter.  Bottom  panel  shows  the  front  view 
of  the  tank  (drawn  to  scale)  with  inside  dimensions  30  cm  x  30  cm  x  15  cm.  This  tank  is  filled 
with  water  to  which  milk  is  added.  The  target  is  suspended  in  the  middle  of  the  tank. 

Fig.  2.  Application  of  the  PDI  system  to  an  aluminum  target  suspended  in  diluted  milk  as 

4 J^*ht 

illustrated  in  Fig.  1;  5  ml  of  milk  was  added.  A  and  B  present  the  images  I||(x,y)  and  I±(x,y) 
convolved  digitally  with  a  two-dimensional  low-pass  spatial  filter.  This  filter  is  roughly  a 
truncated  gaussian,  with  maximal  linear  extent  about  5  pixel  widths;  the  images  of  the  abraded 
patches  are  «  80  x  92  pixels.  The  filter  removes  a  periodic  high  spatial  frequency  artifact 
produced  by  tire  imaging  system.  Image  intensities  (ordinates)  are  expressed  in  the  units  of  the 
8-bit  display,  i.e.,  vary  from  0  to  255  =  28-l.  C  and  D  present  psI(x,y)/2  and  pdl(x,y)+128, 
respectively,  (cf.  Eqs.  1(a)  and  1(b));  an  offset  of  128  was  used,  since  most  pixel  values  of 
pdl(x,y)  are  near  zero,  and  can  be  either  positive  or  negative.  E  and  F  re-present  the  data  of  C 
and  D,  scaled  with  affine  transformations  that  take  advantage  of  the  full  display  range.  The 
transformed  values  are  given  by  ps^y)*^  •  aps[psI(x,y)-pSI(x,y)min]  and  by  pdlfoyltrans  = 
ctpd[pdl(x,y)-pdl(x,y)min]  where  Opg  *  255/[p8I(x,y)max-pSI(x,y)Inin],  and  similarly  for  a^.  In  E, 
psl(x,y)max  and  psl(x,y)min  were  obtained  from  the  disk  region  only;  the  resultant  scale  factor. 
Op,  *  6.4,  is  such  that  the  intensity  variation  of  the  disk  region  occupies  the  full  display  range.  In 
F,  pdl(x,y)max  and  pdl(x,y)mjn  were  obtained  from  the  entire  image,  yielding  apd  «  38.4. 
Abraded  patches,  not  visible  in  Panels  A,  B,  C,  and  D,  are  clearly  visible  in  Panel  F,  but  barely 
visible  in  Panel  E  The  vertical  bands  of  pixels  between  the  arrows  at  the  right  side  of  the  images 
in  A  -  F  were  averaged  over  y  to  generate  numerical  plots  (as  a  function  of  x)  shown  in  panels  A’ 
-F. 
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FIGURE  1 
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The  dose  apposition  of  the  inner  sefmenta  of  the  two  cones  that  combine  to  form  a  double  cone  causes  the  pair 
of  cone  inner  segment!  to  guide  light  as  a  unitary  structure  whose  transverse  sections  are  roughly  elliptical. 
Electron  micrographs  of  the  photoreceptors  of  a  green  sunfish  tLspomie  cyanellus)  retina  provide  evidence  that 
the  refractive  index  in  the  sllipaoid  region  of  the  inner  segments  of  the  double  cones  is  higher  in  the  center  than 
at  the  perimeter.  The  hypothesis  that  the  shape  and  rsfractive-iadex  gradient  could  confer  differential  polar¬ 
ization  sensitivity  on  double  cones  is  examined  with  a  two-dimensional  waveguide  model  of  a  double-cone  inner 
segment.  The  model  has  a  dielectric  constant  that  varies  paraboiicaliy  along  the  narrowest  (x)  dimension,  lead¬ 
ing  to  the  index  profile:  n(x)  -  Hm.I1  -  (x/x«n  ,  when  is  the  peak  value  of  the  index  and  z.  is  a  pa¬ 
rameter  specifying  the  rate  at  which  the  index  decreases  with  increasing  |x|.  A  quantity,  the  polarisation 
contrast,  is  introduced  as  a  measure  at  the  differential  polarisation  sensitivity  of  adjacent  receptors  in  the 
square  mosaic  of  double  cones  in  the  sunfish  retina.  Polarisation  contrast  is  proportional  to  the  relative  dif¬ 
ference  in  power  abeorbed  by  two  double  cones  oriented  with  their  shortest  axse  orthogonal  to  each  other  and 
stimulated  by  a  field  of  uniform  polarisation.  Polarisation  contrast  is  computed  as  a  function  of  wavelength 
for  appropriate  valusa  of  Hm>  and  X*.  For  normally  incident  light  polarised  parallel  to  one  of  the  two  axes  of  the 
double  cones’  cross  sections,  the  polarisation  contrast  is  generally  between  1%  and  6*  for  wavelengths  ranging 
from  #50  to  750  nm.  Over  most  of  those  wavelengths  the  polarisation  contrast  of  the  graded-index-model 
double  cone  is  approximately  five  times  as  large  as  that  of  a  homogeneous-slab  modal  of  the  earns  sise  and  aver¬ 
age  refractive  index.  Additional  benefits  of  a  graded  index,  optical  isolation  Of  adjacent  photoreceptors  and 
antireflection  at  the  photoreceptor  entrance,  are  also  discussed. 


1.  INTRODUCTION 

A.  Vertebrates'  Use  of  Polarisation 
Many  vertebrates  are  capable  of  extracting  information 
from  the  state  of  polarisation  of  visible  light,1"*  but  the 
physical  mechanism  underlying  this  ability  is  not  yst 
established.  In  particular,  the  sensitivity  of  the  green 
aunflsh  l Lepomis  cyanellus)  to  small,  linearly  polarized 
targets  has  been  reported  to  be  a  periodic  function  of  po¬ 
larisation  angle  with  90*  periodicity  and  modulation  depth 
of  aa  mud)  aa  half  a  log  unit*  The  action  spectrum  of 
this  differential  sensitivity  corresponds  to  that  expected 

for  the  pigment  of  the  green  sunfish’s  double  cones  (A - 

621  nm).7  Double  cones*  are  erne  pairs  that  are  fused 
along  the  adjacent  faces  of  their  inner  segments;  they  are 
found  in  the  retinas  of  many  fish,  reptiles,  and  birds*  and 
at  least  one  marsupial10  and  one  monotreme  mammal.* 
In  this  paper  we  examine  the  hypothesis  that  the  shape 
tor,  more  specifically,  the  refractive-index  profile)  of  the 
double  cones  of  sunfish  confers  on  those  cones  a  differen- 
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tial  polarisation  sensitivity.  The  next  subsection  pro¬ 
vides  the  background  for,  and  a  further  elaboration  of, 
this  hypothesis. 

B.  Dodds  Crass  As  Unitary,  Birefringent  Waveguides 
The  index  of  refraction  inside  vertebrate  photoreceptor*  is 
higher  than  the  index  of  the  intracellular  medium  sur¬ 
rounding  them.u  Because  of  this  difference  in  refractive 
indices,  photoreceptors  behave  as  miniature  fiber-optic 
cables,11  guiding  light  from  their  inner  segments  into  (and 
through)  their  outer  segments,  where  phototransduction 
occurs.  Figure  1  shows  a  radial  section  of  the  photorecep¬ 
tor  layer  of  a  green  sunfish  retina  and  illustrates  a  double 
cone  in  longitudinal  section;  the  lower  part  of  the  section 
(dashed  line  marked  b)  is  ~-3  pm  beyond  the  external  lim¬ 
iting  membrane.  Figure  2  shows  a  series  of  transverse 
sections,  taken  at  the  levels  marked  by  the  respective  low¬ 
ercase  letters  in  Fig.  1.  Figures  1  and  2  illustrate  that 
the  two  cones  that  combine  to  form  s  double  cone  are  in 
direct  contact  with  each  other  from  the  external  limiting 
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Pig.  1.  Electron  micrograph  of  the  photoreceptor  layer  of  t  green 
sunfish.  The  sunfish  waa  anesthetized  by  immaraion  in  a  0.2% 
aohitkm  of  tricaine  methanesulfonate  (Sigma).  After  breathing 
had  atopped,  the  heart  waa  perfuaed  with  a  1%  aolution  of  aodium 
nitrite  in  0.05-M  aodium  cacodytate  buffer  to  expel  blood  cella 
from  the  circulatory  lyatem;  thia  waa  followed  fay  perfusion  with  a 
aimilarly  buffered  aolution  of  3%  giutaraldehyde.  The  eyea  were 
removed  horn  the  fiah;  the  lena,  iria,  and  cornea  were  removed; 
and  the  eye  cupa  were  stored  overnight  in  fresh  fixative.  The  eye 
rape  were  rinsed,  dehydrated,  and  stained  in  block  overnight  with 
2%  uranyl  acetate  in  0.05-M  maleate  buffer.  They  were  then 
dehydrated  through  a  graded  aeries  of  alcohol  and  embedded  in 
epon.  Semithin  sections  (1-2  jun  thick)  were  stained  with  Tohii- 
dine  Blue  and  examined  light  microscopically.  Ultrathin  sections 
(0.1  iitn  thick)  were  rat,  mounted  on  formvar  film-coated  (me- bole 
capper  grids,  stained  with  lead  citrate,  and  viewed  in  a  Philips 
300  electron  microecope.  The  plane  of  section  is  perpendicular 
to  the  surface  of  the  retina:  parallel  to  the  long  axis  of  the  photo¬ 
receptors.  Light  entera  from  below  (the  vitreal  direction).  The 
photoreceptor  nuclei  lie  below  the  Held  of  view.  The  dashed  lines 
indicate  the  levels  at  which  the  sections  shown  in  Figs.  2b-2f  were 
rat  The  central  paired  structure  with  the  long  outer  segments, 
os,  is  the  double  cone,  flanked  on  the  right  by  another  double  cone 
and  on  the  left  by  a  single  cone.  The  subceUular  components  are 
identified  in  Fig.  2.  - 


section  are  more  effectively  guided  than  waves  polarized 
parallel  to  the  short  axis,  that  is,  light  polarized  along  the 
long  axis  has  a  greater  ability  to  stimulate  phototrans¬ 
duction.14  Such  geometric  birefringence  could  allow  the 
double  cones  to  act  as  minature,  although  imperfect,  polar¬ 
ization  analyzers:  for  a  given  incident  intensity,  the  frac¬ 
tion  of  power  guided  into  the  photoreceptor  outer  segments 
could  be  a  function  of  the  linear  polarization  state  of  the 
incident  beam. 

Aside  from  depicting  the  shape  of  the  double  cones' 
cross  sections,  Figs.  2  and  3  show  that  the  double  cones  of 
the  sunfish  are  locally  arranged  in  a  roughly  square  mo¬ 
saic, 11  an  arrangement  shared  by  many  teleGals1*-1’  and  at 
least  one  bird.1*  In  this  mosaic  the  long  axes  of  the 
double  cones’  cross  sections  form  the  sides  of  the  squares. 
If  the  hypothesis  that  double  cones  are  geometrically  bi- 
refringent  waveguides  is  correct,  then  the  two  orientation 
classes  of  double  cones  in  the  square-pattern  mosaic  are 
equivalently  two  polarization  classes  of  photoreceptor. 
Bernard  and  Wehner  have  shown  that  polarization  vision 
can  be  understood  in  many  ways  by  analogy  to  color  vi¬ 
sion.1*  In  this  analogy  two  orientation  classes  of  double 
cones  in  a  polarization-sensitive  animal  would  correspond 
to  two  chromatic  classes  of  single  cones  whoee  pigments 
are  maximally  sensitive  to  different  wavelengths  in  a 
color-sensitive  animal.  We  propose  that  the  two  orienta¬ 
tion  classes  of  double  cone  in  the  green  sunfish  are  ideally 
situated  to  function  as  the  input  for  the  computation  of  a 
retinal  polarization-difference  image,  much  as  the  photo¬ 
receptors  of,  for  example,  the  primate  fovea  provide  the 
basis  for  the  computation  of  a  chromatic-difference  im¬ 
age.10  The  significance  of  this  sensory  cell  arrangement 
will  be  discussed  below,  where  it  will  be  clear  that  the 
magnitude  of  the  differential  polarization  sensitivity  of 
the  animal  need  not  be  explained  by  an  equivalent  magni¬ 
tude  of  differential  sensitivity  in  the  individual  receptors. 

To  recapitulate,  in  this  section  we  have  outlined  the 
double-cone- mosaic  hypothesis  of  polarization  vision,  a  hy¬ 
pothesis  that  comprises  two  propositions:  (1)  the  double¬ 
cone  inner  segment  functions  as  a  weak  polarization 
analyzer  by  virtue  of  geometric  birefringence,  and  (2)  the 
mosaic  of  double  cones  serves  as  the  input  to  a  neural  dif¬ 
ferencing  network  that  computes  a  local  polarization  dif¬ 
ference  at  each  point  in  the  animal’s  retina.  Our  focus  in 
this  investigation  is  on  proposition  (1). 


membrane  to  the  bases  of  the  two  outer  segments.  Be¬ 
cause  the  inner  segments  of  the  two  members  of  a  double 
cone  are  in  such  close  apposition,  the  composite  structure 
(i.e.,  the  double-cone  inner  segment)  can  be  expected  to 
act  as  a  single  waveguide.  Evidence  for  energy  coupling 
between  the  two  component  cones  of  a  goldfish  {Carataiut 
auratus)  double  cone  has  been  reported  by  Tobey  et  al.u 

In  many  animals  with  double  cones  the  inner  segments 
of  the  double  cones  have  roughly  elliptical  cross  sections 
in  planes  perpendicular  to  the  guiding  axis.  That  the 
croes  sections  of  the  double  cones  of  the  green  sunfish 
are  roughly  elliptical  can  be  seen  from  Figs.  2  and  3.  We 
hypothesise  that  their  noncircular  cross  sections  and 
material-density  profiles  cause  double  cones  to  exhibit 
geometric  birefringence  such  that  waves  linearly  polar¬ 
ised  parallel  to  the  long  axis  of  the  inner  segment  cross 


2.  GRADED-INDEX  MODEL  OF  THE 
SUNFISH  DOUBLE-CONE  INNER  SEGMENT 

A.  Anatomical  Evidence  That  Double  Cones  are 
Waveguides  with  an  Inhomogeneous  Core 
Figures  1  and  2b-2e  show  that,  as  in  other  vertebrates," 
the  sun  fish's  photoreceptor  inner  segments  have  a  specific 
region  that  is  densely  packed  with  mitochondria.  Photo¬ 
receptor  mitochondria  are  well  known  for  their  essential 
role  in  metabolism  and  are  well  positioned  to  serve  the 
high  energy  demands  of  the  photoreceptor  dark-current 
and  phototransduction  cascade.  However,  the  concentra¬ 
tion  of  mitochondria  in  vertebrate  photoreceptors  also 
serves  to  increase  the  photoreceptors’  refractive  index,  en¬ 
hancing  the  inner  segments’  ability  to  guide  light.11  We 
now  concentrate  our  attention  on  the  mitochondria  be¬ 
cause  of  their  role  in  the  establishment  of  waveguiding  in 
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Fig.  2.  Six  electron  micrographs  at  different  levels  through  the  same  tetrad  of  cones,  all  at  the  same  magnification  and  orientation.  The 
preparation  of  the  tiasue  was  the  ax=ne  as  that  for  Fig.  1.  The  plane  of  section  is  the  same  in  all  panels:  parallel  to  the  retinal  surface 
and  perpendicular  to  the  long  axis  of  the  photoreceptors;  the  approximate  levels  of  panels  b-f  are  shown  in  Fig.  1.  a:  This  section  passes 
through  the  photoreceptor  nuclei  (at  a  level  below  the  Held  shown  in  Fig.  1).  The  rod  nuclei,  m,  are  finely  granular;  the  cone  nuclei,  cn, 
an  more  mottled.  The  nine  cones  that  form  one  square  of  the  double-cone  mosaic  are  indicated  by  the  numbers  1— St  and  some  or  all  of 
these  cooes  are  numbered  in  other  panels  as  well,  b:  This  section  passes  through  the  base  of  the  cone  inner  segments;  the  paraboloid,  p, 
and  the  mitochondrial  cluster,  m,  are  indicated  in  the  four  double  cones.  They  are  surrounded  by  the  fine  processes  of  the  Mttller  cells,  me. 
and  the  rod  myoids,  rm,  the  attenuated  processes  that  connect  the  rod  nuclei  to  their  inner  segments.  The  arrowheads  indicate  the  dark 
intercellular  junction  between  the  two  members  of  the  double  cones,  c;  At  this  more-distal  level  the  cone  inner  segments  occupy  essen¬ 
tially  the  entire  cross  section  of  the  retina,  and  the  cross  section  of  the  interior  of  the  double  cone  is  nearly  filled  by  the  mitochondrial 
duster.  The  electron  density  of  the  mitochondrial  cluster  is  greater  than  in  b,  and  the  gradients  of  electron  density  across  the  double 
cones  are  evident,  d:  At  this  level  of  section  the  cones  fill  most  of  the  retinal  cross  section,  but  pigment  granules,  pg,  from  the  pig¬ 
mented  epithelium  begin  to  intrude  from  above.  Within  each  double  cone  the  mitochondrial  cluster  occupies  a  smaller  fraction  of 
the  cross  section  than  in  c,  reflecting  the  taper  evident  in  Fig.  1.  The  gradient  in  electron  density  is  still  evident,  e:  The  mitochon¬ 
drial  dusters  have  tapered  still  more,  the  pigment  granules  are  more  numerous,  and  the  extracellular  space  around  the  cones  has  in¬ 
creased  relative  to  that  shown  in  d.  f:  The  cone  outer  segments,  os,  occupy  a  much  smaller  fraction  of  the  retinal  cross  section;  the 
pigmented  processes,  a  much  larger  one. 
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Fig.  3.  Semithin  sections  prepared  as  described  in  the  caption 
of  Fig.  1  and  viewed  with  a  light  microscope.  This  figure  demon¬ 
strates  the  regularity  in  the  photoreceptor  mosaic  of  the  green 
sunfish.  The  regularity  seen  here  is  representative  of  the  mosaic 
across  the  entire  retina  (see  Ref.  15) 


photoreceptors.  Figures  1  and  2  illustrate  three  features 
that  we  believe  are  of  fundamental  importance  for  under¬ 
standing  the  mitochondria’s  influence  on  the  waveguiding 
behavior  of  double-cone  inner  segments. 

Mitochondria  concentrated  near  the  double-cone  center. 
First,  the  distribution  of  mitochondria  in  each  of  the  two 
cones  that  unite  to  form  a  double  cone  is  concentrated  near 
the  common  boundary  of  the  two  (i.e.,  along  the  fusion 
boundary).  This  asymmetrical  packing  of  mitochondria 
in  the  individual  cones  that  combine  to  form  a  double  cone 
stands  in  sharp  contrast  to  the  distribution  of  mitochon¬ 
dria  in  the  single  cones  seen  at  the  center  of  each  square 
of  the  mosaic  in  Fig.  2.  In  the  single  cones  the  mitochon¬ 
drial  cluster  is  centered  in  the  cytoplasm  at  a  more-or-less- 
uniform  distance  from  the  cytoplasmic  membrane. 

Elliptical  annulus  of  cytoplasm  surrounding  mito¬ 
chondria.  A  second  feature  of  the  distribution  of  mito¬ 
chondria  that  we  believe  to  be  a  determinant  of  the 
waveguiding  behavior  of  the  double-cone  inner  segments  is 
that  the  region  containing  mitochondria  is  surrounded  by 
an  annulus  of  cytoplasm  devoid  of  mitochondria.  Thus 
the  mitochondrial  region  is  contained  within  a  reasonably 
well-defined  boundary  that  has  the  same  shape  as  the 
outer  envelope  of  the  double  cone  but  is  clearly  interior  to 
that  outer  boundary.  Interestingly,  there  is  one  perspec¬ 
tive  from  which  the  distribution  of  mitochondria  in  the 
double  cone  resembles  that  in  the  single  cone:  each  has 
an  annulus  of  cytoplasm  circumscribing  the  mitochon¬ 
drial  region.  In  the  case  of  the  double  cone,  however,  the 
annulus  of  cytoplasm  follows  the  perimeter  of  the  compos¬ 
ite  structure  rather  than  the  perimeter  of  the  individual 
member  cones.  That  is,  the  distribution  of  mitochondria 
within  a  single  cone  resembles  the  distribution  of  mito¬ 
chondria  in  a  double  cone  as  a  composite  structure  and  not 
the  distribution  of  mitochondria  in  an  individual  member 


inner  segment,  while  comprising  two  distinct  photorecep¬ 
tors,  appears  from  the  optical  perspective  to  have  a  single 
ellipsoid.  That  is,  the  mitochondrial  distribution  in  the 
double-cone  inner  segment  forms  a  single  functional  body 
that  we  may  call  “the  double-cone  ellipsoid. "  We  consider 
the  existence  of  the  unitary  double-cone  ellipsoid  to  be  an 
adaptation  that  functions  to  unify  the  member  inner  seg¬ 
ments  of  a  double  cone  into  a  single  waveguide. 

Greater  electron  density  near  the  centers  of  ellipsoid 
cross  sections.  The  third  characteristic  feature  of  the 
mitochondrial  distribution  in  double-cone  inner  segments 
is  that  the  mitochondria  nearest  the  center  of  transverse 
sections  of  the  ellipsoid  are  more  electron  dense  than  the 
mitochondria  near  the  ellipsoid's  perimeter.  Figures  2c- 
2e  show  a  clear  gradient  in  electron  density  across  the  mi¬ 
tochondrial  cluster  in  each  double  cone.  This  gradient  in 
electron  density  cannot  represent  a  differential  packing  of 
mitochondria  along  the  *  axis  (into  and  out  of  the  micro¬ 
graphs),  because  these  sections  are  only  0.1  pm  thick, 
much  less  than  the  thickness  of  a  typical  mitochondrion. 
Thus  the  gradient  is  in  the  form  of  a  differentiation  among 
mitochondria:  the  mitochondria  toward  the  center  of  the 
cell  are  more  electron  dense  than  those  at  the  periphery. 

It  seems  unlikely  to  us  that  the  observed  differential 
electron  density  of  mitochondria  is  an  artifact  of  the  pene¬ 
tration  of  the  fixatives  and  stains  used  to  prepare  the  tis¬ 
sue  (as  described  in  the  caption  of  Fig.  1).  In  the  case  of 
a  penetration  artifact,  we  would  expect  the  staining  to  be 
a  simple  decreasing  function  of  the  distance  to  the  nearest 
point  on  the  plasma  membrane.  However,  the  micro¬ 
graphs  show  that  this  is  not  the  case:  the  electron  den¬ 
sity  of  the  mitochondria  is  an  increasing  function  of  the 
mitochondrion’s  distance  from  the  cell  membrane.  More¬ 
over,  at  the  scleral  end  (near  the  outer  segments)  the  mito¬ 
chondria  are  darker  than  at  the  vitreal  end  (near  the 
nucleus)  even  in  the  regions  immediately  adjacent  to  the 
cytoplasmic  envelope  (see  Fig.  1).  Thus  whatever  process 
determines  the  observed  electron  density  of  individual 
mitochondria  apparently  operates  differentially  in  all 
three  dimensions.  The  variation  in  mitochondria  as  a 
function  of  position  in  the  inner  segment  is  not  unique  to 
sunfish.  It  is  well  established  that  photoreceptor  mito¬ 
chondria  are  often  modified,  relative  to  the  mitochondria 
in  other  cells,  in  ways  that  affect  the  optical  properties  of 
the  photoreceptor,  and  these  modifications  are  more 
prominent  around  the  perimeter  and  the  scleral  end  of  the 
ellipsoid.”*" 

The  differential  electron  density  of  the  mitochondrial 
regions  in  the  sections  shown  in  Fig.  2  must  represent  an 
increase  in  the  concentration  of  either  an  electron-dense 
substance  or  a  substance  for  which  the  fixative  or  stains 
have  an  affinity.  The  osmium  tetroxide  used  in  these 
preparations  binds  principally  to  lipids,  the  uranyl  acetate 
binds  to  nucleic  acids,  and  the  lead  citrate  binds  to  both 
osmium  and  uranium  atoms.  If  the  variation  in  electron 
density  arises  from  a  differential  affinity  for  the  stain,  it 
will  still  indicate  a  correlated  variation  in  the  concentre* 


of  a  double-cone  pair. 

In  the  anatomical  literature  the  mitochondrial  region  of 
the  inner  segment  of  single  photoreceptors  is  known  as 
the  ellipsoid.  This  name  refers  to  the  structure’s  three- 
dimensional  shape  in  single  cones.”  7ft  can  summarize 
our  first  two  observations  by  saying  that  the  double-cone 


tion  of  lipid  or  nucleic  acids  in  the  unstained  tissue.  An 
increase  in  the  concentration  of  nucleic  acids  or  lipids  re¬ 
sults  in  an  increase  in  refractive  index.” 

In  summary,  analysis  of  electron  micrographs  sup¬ 
ports  the  general  hypothesis  that  the  refractive  index  of 
the  ellipsoid  region  of  sunfish  double  cones  is  highest  in 
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Dm  canter  of  the  croaa  auction  and  daclinaa  toward  the 
perimeter. 

B.  Mathematical  Analyais  of  the  Double- Cone  Inner 
Segment1!  Waveguiding  Properties 

While  our  anatomical  observations  support  the  proposition 
that  the  double-cone  inner  segment  is  a  unitary  graded- 
index  waveguide  with  a  noncircular  cross  section,  alone 
they  do  not  provide  sufficient  quantitative  constraints  for 
analysing  the  properties  of  the  guide.  In  this  section  we 
synthesise  our  anatomical  observations  with  additional 
information  and  hypotheses  into  a  mathematical  model 
suitable  for  examining  the  waveguiding  behavior  of  the 
double  cone. 

From  the  outset  we  acknowledge  that  we  have  made  sev¬ 
eral  simplifying  assumptions  in  constructing  the  model. 
One  of  the  crudest  of  these  is  that  our  model  neglects 
variation  in  properties  along  the  longitudinal  (z)  axis.  We 
justify  our  neglect  of  z-axis  variation  at  this  $ime  by  the 
need  for  a  relatively  simple,  tractable  model.  Nonetheless, 
the  anatomical  data  that  we  have  presented  show  clear 
and  systematic  variation  along  the  z  axis  that  may  play  an 
important  role  in  the  waveguiding  behavior  of  the  double 
cone.  We  take  up  this  matter  again  in  Section  4. 

In  addition,  in  order  to  investigate  quantitatively  the 
waveguiding  properties  of  the  double  cones,  we  have  made 
specific  assumptions  about  the  shape  and  height  of  the 
refractive-index  profile.  Constraints  on  these  assump¬ 
tions  arise  from  measurements  of  the  size  of  the  double 
cones,  from  previous  estimates  of  the  average  refractive 
index  of  the  ellipsoids  of  photoreceptors,  and  from  the 
maximum  value  that  the  refractive  index  can  attain  in  a 
living  cell. 

Size  of  the  waveguide  cross  section.  From  the  electron 
micrographs  of  Figs.  1  and  2  it  is  clear  that  the  size  of  the 
double-cone  cross  section  varies  somewhat  with  position 
along  the  z  axis.  Nonetheless,  as  noted  above,  we  have 
chosen  to  neglect  z-axis  variations  and  have  selected  a 
single  (x-y)  croaa  section  of  a  double  cone  as  representa¬ 
tive.  In  order  to  generate  a  conservative  estimate  of  the 
polarization  properties  of  the  double  cone,  we  baaed  our 
model  on  the  largest  cross  sections  that  have  been  mea¬ 
sured  in  living  double  cones.  The  dehydration  that  is  a 
necessary  step  in  the  preparation  of  the  retina  for  viewing 
the  tissue  with  an  electron  microscope  causes  the  cells  to 
shrink  somewhat.  Because  of  this  shrinkage,  it  is  better' 
to  estimate  the  size  of  the  cones  by  viewing  them  with  a 
light  microscope  under  more-natural  conditions.  Light- 
microscopic  observations  of  fixed  and  stained  retinas  or  of 
fresh  cones  suspended  in  fish  Ringer’s  solution  indicate 
that  a  resonable  upper  limit  to  the  largest  cross  sections  of 
sunfish  double-cone  inner  segments  is  7  pm  x  14  p.m,15 

Average  refractive  index.  Using  refractometry  and  dif¬ 
ferential  interference  contrast  microscopy,  Sidman  studied 
the  photoreceptors  of  many  animals.11  His  measurements 
showed  that  the  ellipsoid  region  of  cone  inner  segments 
have  an  average  refractive  index  of  approximately  1.3ft 
Although  he  did  not  investigate  green  sunfish  in  his  study, 
he  did  examine  photoreceptors  of  a  wide  array  of  verte¬ 
brates  including  frog,  rat,  chicken,  and  monkey.  Initially 
we  adopt  the  value  1.39  as  the  average  refractive  index 
across  the  double-cone  inner  segment. 

Refractive-index  profile.  On  the  basis  of  the  anatomi¬ 


cal  evidence  presented  above,  we  conclude  that  the  refrac¬ 
tive  index  of  the  double-cone  inner  segment  is  highest  in 
the  center  of  the  structure  and  declines  monotonically  to 
the  perimeter.  However,  the  anatomical  data  cannot  be 
converted  into  a  specific  formula  for  the  refractive-index 
profile,  because  there  is  no  established  relationship  be¬ 
tween  the  electron  density  in  a  fixed  and  stained  retina 
and  the  electron  density  in  living  photoreceptors.  Fig¬ 
ure  4  illustrates  how  our  choice  of  refractive-index  profile 
is  motivated  by,  and  consistent  with,  our  anatomical  obser¬ 
vations.  Figure  4A  is  a  magnified  view  of  cones  8  and  9 
from  Fig.  2c.  Figure  4B  is  a  contour  plot  of  the  digitally 
filtered  density  of  the  negative  used  to  generate  Fig.  4A 
(see  Fig.  4  caption  for  details).  We  embed  the  contour  plot 
in  a  slab  waveguide  to  indicate  the  coordinate  system  used 
throughout  this  paper.  Figure  4C  shows  three  lines  of  the 
digitally  filtered  data  taken  from  the  locations  indicated  by 
the  lines  in  Fig.  4A.  The  mesh  plot  in  Fig.  4D  provides  a 
different  view  of  the  data  presented  in  the  contour  plot  of 
Fig.  4B.  Figures  4C  and  4D  include  a  curve  that  matches 
the  index  profile  that  we  will  be  using  in  our  model.  We 
began  by  approximating  the  variation  of  the  dielectric 
constant  as  a  two-dimensional  parabolic  gradient: 

e(x,y)  -  £nu[l  -  (x/x0)*  ~  (y/yo)*l  ■  (1) 

The  dielectric  constant  is  proportional  to  the  square  of  the 
refractive  index,  so  Eq.  (1)  leads  to 

n(x,y)  -  rw[1  -  (x/x0)*  -  (y/yo)*]l/*.  (2) 

The  index  variation  of  Eq.  (2)  is  consistent  with  our  ana¬ 
tomical  data  and  has  some  precedent  in  the  literature.*7-** 
Equations  (1)  and  (2)  are  assumed  valid  for  the  intracellu¬ 
lar  region  of  the  photoreceptors,  which  can  be  approxi¬ 
mated  by  the  region  in  which  [x/(d/2)]*  +  (y/d)’  s  1, 
where  d  is  the  length  of  the  minor  axis  of  the  model  double 
cone’s  cross  section.  The  constants  e**,  and  n^  are  the 
values  of  the  dielectric  constant  and  the  refractive  index 
at  the  center  of  the  guide.  The  constants  xo  and  yo  deter¬ 
mine  bow  rapidly  the  dielectric  constant  and  the  refractive 
index  drop  with  distance  away  from  the  origin  along  the  x 
and  the  y  axes,  respectively.  The  parameters  x0  and  y0 
can  also  be  expressed  in  terms  of  the  above  parameters  d 
and  a,*,-.  y0  “  2xo  “  d[Vl  -  (n£//n1M,)l)'1'J,  where  nc, 
is  the  index  of  the  extracellular  fluid.  With  d  »  7  pm, 
nd  -  1.33,  and  with  Sidman’s  value  of  1.39  for  the  average 
value  of  the  refractive  index  across  the  inner  segment  at 
the  ellipsoid  region,11  the  values  of  the  parameters  nM„ 
x0,  and  y0  are  determined  and  readily  computed:  - 

1.42,  Xo  -  10.0  pm  andy0  -  20.0  ^im. 

In  constructing  the  model  we  made  one  further  approxi¬ 
mation.  We  chose  to  neglect  the  variation  in  refractive 
index  along  they  axis,  so  that  the  assumed  index  profile  of 
our  model  is  one  dimensional: 

nix)  -  n*tt[l  -  (x/xo)*],/*  (3) 

The  neglect  of  y-axis  variation  considerably  simplifies  the 
calculations  and  is  not  unreasonable  given  that  the  varia¬ 
tion  along  the  x  axis  dominates  the  behsvior  of  the  guide.** 

Bound-mode  solutions  of  the  model  waveguide.  To 
solve  for  the  modes  of  our  model  double  cone,  we  used  the 
infinite  parabolic  profile.  This  model,  where  Eq.  (3)  is 
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^*8-  between  the  anatomy  and  the  mathematical  model.  A,  A  reproduction  of  one  of  the  double  cones  (8-9)  from  the 

section  shown  in  Fig.  2c.  The  scale  bar  at  the  upper  left  of  the  cell  is  620  nm  long;  however,  it  must  be  remembered  that  this  tissue  has 
been  shrunk  somewhat  by  its  preparation.  The  other  lines  refer  to  C,  below.  B,  The  negative  that  generated  A  was  photographed  against 
a  backlit  piece  of  white  acrylic  by  a  Hamamatsu  CCD  camera  whose  output  went  to  an  Imaging  Technologies  Series  161  image-analysis 
system.  After  we  had  calibrated  the  camera/image-analysis  system  to  ensure  linearity  of  the  output  over  intensity  and  space,  we  used  the 
series  161  to  digitise  the  image  of  the  negative  as  acquired  by  the  CCD  camera.  The  digital  image  data  was  then  loaded  into  the  matlab 
computer  program  (version  4.0,  the  MatWforks,  Inc.).  In  matlab  the  electron-density  data  was  convolved  with  a  Gaussian  filter  that  ap¬ 
propriately  averaged  the  data  to  convert  the  density  profile  to  an  approximation  of  the  index  profile.  The  filter  fell  in  value  to  10%  of  its 
maximum  at  820  nm  from  the  canter  and  was  normalised  by  its  integral  over  all  x  and  y  being  set  to  1.  The  filtered  image  data  were  then 
merited  so  that  only  the  region  within  the  cytoplasm  waa  displayed,  and  a  contour  plot  was  generated.  The  contour  plot  is  viewed  at  an 
angle  and  truncated  to  fit  into  a  slab  indicating  the  coordinate  system  used  in  the  model.  Note  that  the  contour  lines  follow  the  general 
outline  of  the  compoeite  structure;  that  optically  the  boundary  between  the  ellipeoid  of  cell  8  and  the  ellipsoid  of  cell  9  is  essentially 
nonexistent.  C,  Three  lines  from  the  filtered  image,  data  were  scaled  to  suggest  the  form  of  the  refractive-index  gradient.  We  accom- 
plished  scaling  by  first  setting  the  index  of  the  extracellular  regions  to  1.33  and  then  setting  the  average  density  across  the  x  axis  near  the 
center  of  the  cell  to  1.39.  The  three  lines  that  we  scanned  to  generate  this  plot  are  indicated  in  A.  The  dotted  curve  came  from  the  top 
Une  in  A,  and  the  dotted-dashed  curve  came  from  the  crater  line  of  A  (only  the  ends  are  drawn,  to  as  not  to  obscure  the  details  of  the 
figure).  The  dashed  curve  is  the  bottom  line  in  A.  The  solid  curve  is  a  least-squares  fit  of  the  data  from  the  top  line  in  A  (the  dotted 
curve  of  C).  The  solid  curve  satisfies  Eq.  (3)  with  —  1.42  and  xo m  8  nm  (leading  to  d  —  6.6  nm).  We  stress  that  the  numbers  along 
the  ordinate  are  crude  approximations  that  do  not  derive  from  direct  measurements  of  the  refractive  index.  D,  Mesh  plot  of  the  data 
displayed  in  B.  The  solid  curve  used  to  fit  the  data  in  C  is  drawn  at  10  times  the  thickness  of  the  grid  lines  of  the  mesh. 


valid  for  all  x,  has  the  virtue  that  the  bound-mode  solu¬ 
tions  of  the  vector  wave  equation  derived  from  Max  trail's 
equations  and  Eq.  (3)  are  in  analytical  closed  form  for  inci¬ 
dent  radiation  polarized  along  they  axis.  Moreover,  with 


one  more  approximation  (presented  below),  there  exists  an 
analytical  solution  for  the  bound  modes  excited  by  radia¬ 
tion  polarized  along  the  x  axis,  as  well.  The  solution  of 
the  equations  for  the  guided  modes  excited  by  light  polar- 
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food  along  the  y  uii  yield*  the  TE  modes,"  which  have 
electric  field*  of  the  form 

£,„(«)-  ATi.«.(«)exp(-5a/2),  (4) 

where  6  »  (Ao/xo)1'3*.  *«  is  the  free-space  wave  num¬ 
ber,  Art.  «  a  parameter  used  to  normalize  modes  such 
that  each  carries  unit  power  per  unit  length  of  the  y  axis, 
and  Hm(6)  is  the  Hermite  polynomial  of  order  m  (m  - 
0, 1, 2, . . .  ).31 

Since  we  will  be  assuming  a  monochromatic  incident 
field  with  angular  frequency  <u,  and  since  we  are  ignoring 
the  variations  in  size  and  index  along  the  z  axis,  all  fields 
of  interest  will  be  waves  whose  z  and  t  dependence  is  ex¬ 
pressed  in  the  factor  exp[t(<o*  -  /9te„*)]-  It  can  be 
shown*1  that  Sts.  must  satisfy 

a  *  /■  ,i  (2m  +  Dio^nu  /p, 

0TI«  •  (*0«««) - (5) 

*0 

In  the  second-order  differential  equation  that  must  be  sat¬ 
isfied  for  an  electromagnetic  wave  propagating  along  the 
z  axis  and  polarized  along  the  x  axis  (the  TM  modes),  we 
can  discard  the  terms  of  order  x*/x*-,  those  terms  would 
change  the  coefficient  of  the  electric  field  in  this  equation 
by  at  moat  0.0023%  for  —  1.42  and  x0  -  10.0  pm. 
The  solutions  to  the  equation  without  those  terms  are  of 
the  form 


perimeter  of  our  model  double  cone.  Figure  5  shows  the 
form  of  all  the  even-symmetric  modes  that  meet  this  crite¬ 
rion  when  nMI  1.42  and  *o  ■*  10.0  pm.  The  value  of  Xo 
determined  in  the  caption  for  Fig.  4  is  smaller  than 
10.0  pm,  but,  as  discussed  above,  we  begin  our  analysis  with 
a  limiting  case.  In  Fig.  5  we  show  only  even-symmetric 
modes,  because  given  the  incident  radiation  that  we  as¬ 
sume  below,  these  are  the  only  modes  that  will  be  excited. 
Note  the  magnitudes  of  the  fields  where  x  -  d/2. 

Calculation  of  the  power  extracted  from  different  input 
polarizations.  We  now  examine  the  power  extracted  by 
our  model  photoreceptor  when  it  is  excited  by  a  uniformly 
polarized,  normally  incident  plane  wave-33  We  assumed 
that,  before  the  entrance  of  the  waveguide,  the  index  is 
homogeneous  and  equal  to  that  of  the  extracellular  fluid. 
The  change  in  refractive  index  at  the  interface  leads  to  a 
reflection,  which  we  modeled  by  calculating  the  Fresnel 
transmission  coefficient  at  each  point  along  the  cross  sec¬ 
tion.34  We  then  multiplied  the  amplitude  of  the  incident 
plane  wave  by  the  resulting  function  of  x.  Use  of  the 
Fresnel  coefficients  is,  of  course,  another  approximation, 
but  its  simplicity  makes  its  usage  common  M  The  guided 
field  can  then  be  determined  by  expansion  of  the  transmit¬ 
ted  field  into  the  modes  described  by  Eqs.  (4),  (6),  and  (7), 
since  the  modes  form  an  orthogonal  set  over  the  infinite 
cross  section.  That  is,  the  electric  field  of  any  guided 
wave  can  be  expressed  as  the  sum  of  these  modal  solutions: 


E,,m(p) 


At«, 

n(x) 


ff„(p)exp 


(6) 


where  p  ■  x[(*0nw/x0)3  +  4/x04]u*  and  Atm.  i>  a  normal¬ 
isation  parameter  analogous  to  Atc.. 

In  the  TM  case  there  is  a  small  longitudinal  component 
of  the  electric  field  that  can  be  found  from  the  transverse 
component  with  the  equation 


„  ,  ,  -i  fd(E,..(x)]  2xE,.*(x)l 

01 


where  /3tm.  ia  determined  by 


0TM.*  -  (A.*™,)3  “  ~5  "  (2m  +  l)[(~fi)i  +  £]" 


(8) 


As  indicated  above,  the  fields  described  by  Eqs.  (4M8)  are 
exact  (within  the  approximation  for  the  TM  modes)  only  if 
Eq.  (3)  is  valid  for  all  x.  However,  they  represent  a  good 
approximation  to  the  fields  bounded  by  a  guide  described 
by  Eq.  (3)  valid  for  the  region  inside  the  double  cone  if  all 
the  field  components  [e.g.,  from  Eqs.  (4),  (6),  and  (7)]  die 
off  quickly  before  reaching  the  region  over  which  Eq.  (3)  is 
no  longer  valid.  In  order  to  use  the  solutions  represented 
by  Eqs.  (4)-(8),  we  adopted  the  criterion  that  a  mode  be 
considered  bound  only  if  its  effective  index  was  greater 
than  1.33.w  A  physical  intuition  of  the  significance  of 
this  restriction  can  be  obtained  from  the  WKBJ  approxi¬ 
mation.  In  that  analysis,  bound  waves  are  reflected  back 
toward  the  center  of  the  guide  at  the  point  where  the  re¬ 
fractive  index  is  equal  to  the  effective  index.  Our  restric¬ 
tion  is  equivalent  to  one  in  which  a  mode  is  considered 
guided  only  if  the  so-called  walls  of  the  guide  are  inside  the 


Eguidtd  “  2)  Re{pTM.E,..(x)exp[t(wt  -  0tm.z  +  4>tm.)]}* 

m 

+  R*{pn_  E*»(x)exp[i(«t  -  0-ti.z  +  4>-rO]}y 

+  Re{pTM.£fc»(*) 

x  exp[i(a>f  —  0tm.z  +  d>Tum  ~  w/2)]}z ,  (9) 

where  Re  means  the  real  part  of,  x,  j>,  and  f  are  unit  vec¬ 
tors  along  the  x,  y,  and  z  axes,  respectively,  the  p’s  are  am¬ 
plitude  coefficients,  the  <£’s  are  relative  phases  of  the 
modal  fields,  and  all  other  terms  are  as  defined  above. 

The  polarization  of  any  normally  incident  wave  can  be 
resolved  into  two  waves  polarized  along  the  x  and  y  axes, 
so  it  is  sufficient  to  analyze  the  propagation  of  the  fields 
resulting  from  incident  waves  at  these  two  polarizations 
only.  Furthermore,  since  a  normally  incident  plane  wave 
has  even  symmetry  around  the  plane  defined  by  x  —  0,  we 
can  determine  by  inspection  that  the  wave  will  excite  only 
even-indexed  modes. 

We  thus  obtain  five  bound  even  modes  for  each  of  the 
two  incident  polarizations.  The  relative  contribution 
of  each  of  the  modes  was  determined  by  numerical  in¬ 
tegration  with  software  written  by  us  in  C,  and  also 
with  Mathematics  (version  1.2,  Wolfram  Research,  Inc.). 
These  values  can  be  found  in  Table  1.  In  order  to  deter¬ 
mine  what  these  numbers  mean  in  terms  of  excitation  of 
the  double  cones,  we  calculated  the  fraction  of  propagated 
energy  that  would  be  available  for  driving  phototransduc¬ 
tion  (i.e.,  chromophore  isomerization).  We  will  continue 
to  use  the  modes  determined  from  the  lossless  parabolic 
profile  model,  although  we  recognize  that  for  the  outer 
segments,  basically  a  pair  of  tapering  circular  cylinders, 
this  approximation  is  somewhat  coarse.  The  analysis  will 
at  least  provide  a  reasonable  approximation  to  the  fields 
present  at  the  entrance  to  the  outer  segments,  because 
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Fig.  5.  Form  of  tranevarae  compononU  of  electric  Add*  that 
propagate  to  the  panbolic-indaa-profile  modal.  For  the  values  of 
Om  and  «f  comepooding  to  the  analysis  in  Table  1,  tbooo  an  the 
ontjr  modao  that  an  both  bound  (have  affactiva  indicaa  greater 
than  1.33)  and  capable  of  excitation  by  a  normally  incident  uni- 
form  plana  wave.  The  electric  field  excited  by  ay -polarized  inci¬ 
dent  beam  ie  conatrained  more  closely  to  the  cantor  of  the  guide 
(a  -  0),  but  at  thia  aeale  the  field  diatributiooa  an  indiatinguiah- 
able.  Depicted  ban  an  the  modee  for  m  -  0,  2,  4,  8,  and  8  in 
bp.  (4)  *t  ttq.  Bean  when  m  —  8,  the  field  die*  out  quickly  be¬ 
yond  the  walla  of  the  guide  (when  x  -  3.5  pm). 


Table  I.  Coefficient!  Obtained  from  Analysis  of  the  Graded-Index  Model* 


Coefficient 

Mode  Number 

0 

2 

4 

6 

8 

An  (V/m) 

13420 

4778 

694.7 

63.89 

4.302 

An  (V/m) 

18980 

6663 

964.7 

86.60 

6.736 

0n  (pm*1) 

14.3347 

14.1329 

13.9283 

18.7206 

13.6097 

0tm  (pm'1) 

14.3343 

14.1826 

13.9279 

13.7202 

13.6092 

Pn 

0.0028 

0.0020 

0.0017 

0.00167 

0.0016 

Pvt 

0.0028  ’ 

0.0020 

0.0017 

0.00167 

0.0016 

a«i  (pm’) 

38.48 

19.62 

14.86 

12.62 

10.63 

o«x  (pm*) 

38.39 

19.30 

14.64 

12.12 

10.16 

*Thie  table  rcprceente  the  coefficient*  for  a  waveguide  whose  iadea-of .refraction  profile  ie  eatirted  by  Iq.  (3),  with  x»  — 100  pm  and  -  1.42.  All  the 
numbers  in  tbia  table  ccwaepond  to  the  modee  eeeociatad  with  aa  Incident  field  wheee  wavelength  la  620  am  la  free  epees.  The  A'e  are  coeffldente  that  are 
eatabtWmd  to  ertbenermaUee  the  modee;  with  tboee  ooefficieate  each  mode  carriee  unit  power  per  unit  length  of  they  exie.  The  $'*  ere  the  propagation 
constant!,  or  longitudinal  ware  munbere.  The  p’s  a<a  avpeneie"  coefficients  that  are  obtained  when  e  normally  incident  uniform  plane  wave  la  expanded  Into 
tbo  bound  orthenormal  modes;  the  TE  and  TM  modao  an  cached  when  tho  incMaot  ware  ie  poiarited  along  tbo  y  exie  end  the  x  axie,  respectively.  The  a  * 
an  defined  bgr  Eqa.  (10)  and  (ID,  with  fU-1  W/m*.  gubtrectlng  the  corresponding  a’e  and  dividing  by  the  average  of  tbo  power  received  from  ooch  of  the 
tocidei*  polariaatiene,  one  obtain  Cp  ae  defined  by  Iq.  (12).  For  this  eaae  the  value  ie  1.6*. 
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tho  actual  fidds  should  bo  similar  to  our  model  fields  in 
the  ellipsoids,  and  the  ellipsoids  directly  abut  the  bases 
of  the  outer  segments.  As  can  be  seen  from  Fig.  1, 
the  perimeters  of  the  outer  segment  and  the  ellipsoid  are 
contiguous. 

The  abeorption  of  light  by  a  photopigment  is  propor¬ 
tional  to  the  square  of  the  scalar  product  of  the  electric- 
field  vector  and  the  transition-moment  vector  of  the 
absorbing  chromophore."  For  the  TE  modes  we  need  con¬ 
cern  ourselves  only  with  the  square  of  the  electric  field, 
because  the  pigments  are  free  to  rotate  in  the  disks  of  the 
photoreceptor  outer  segments,*1"  and  thus,  on  average  over 
time,  all  orientations  of  the  chromophore  in  the  x-y  planes 
will  be  represented  equally.  However,  in  the  case  of  the 
TM  modes  there  is  a  small  x  component  of  the  E  field,  so 
the  dichroism  of  the  pigments  must  be  taken  into  account. 
A  conservatively  small  estimate  of  the  axial  dichroism 
of  the  pigments  based  on  an  analysis  of  the  dichroism  of 
photoreceptor  outer  segments  stimulated  from  their  sides 
is  3 : 1."  With  this  information  and  with  a  measure  of  the 
absorbing  ability  of  the  photopigment  we  calculated  the 
power  absorbed  by  a  double  cone.  As  a  measure  of  the  ab¬ 
sorbing  ability  of  the  pigment,  we  used  the  specific  optical 
density  of  the  long-wavelength-sensitive  cones  of  the  gold¬ 
fish  (0.0128  ixm~l  O.D.  units  at  its  Am.  of 625  nm).40  With 
that  value  for  the  specific  optical  density,  we  approximated 
the  rate  at  which  energy  is  absorbed  along  the  z  axis  fay 
presuming  that  the  x  dependence  of  the  rate  of  absorption 
is  carried  only  by  the  dependence  on  x  of  the  field  strength 
[i.e.,  that  the  power  decays  as  exp(- Dz),  where  D  is  a  con¬ 
stant  equal  to  ln(10)  times  the  specific  optical  density  of 
the  pigment].  The  conductivity  of  the  medium  is  thus  D/it, 
where?)  is  the  characteristic  impedance  of  the  volume.  In 
the  determination  of  ij,  we  used  an  index  of  refraction  of 
1.385,  the  approximate  index  of  the  cone  outer  segment.11 
Integrating  the  conductivity  multiplied  fay  the  square  of 
the  electric  field  at  each  point  in  the  volume  gives  the  ab¬ 
sorbed  power.  Division  of  the  absorbed  power  fay  the  inci¬ 
dent  power  density  results  in  an  effective  collecting  area 
for  the  double  cone.  For  TE  modes  (which  are  excited  by 
incident  waves  polarised  along  the  y  axis),  we  obtain 

°"1 “  }  J^UW*>l*d*  •  (10) 

For  TM  modes  (which  are  excited  fay  light  polarized  along 
the  x  axis)  we  obtain 

(id 

where  the  division  by  3  of  the  energy  in  the  longitudinal 
field  component  represents  the  effect  of  the  photopig¬ 
ment's  dichroism.  In  Eqi.  (10)  and  (11),  Pm  is  the  incident 
power  flux  density,  /  is  the  length  of  the  outer  segment, 
and  d  is  the  thickness  of  the  guide- the  distance  between 
which  n  *  1.33— along  the  x  axis.  The  factor  of  2d  outside 
the  integrands  sppears  as  a  consequence  of  the  integration 
from  -d  to  d  along  the  y  axis.  Initially  d  was  fixed  at 
7  Mm,  but  in  later  analyses  we  allowed  it  to  vary  in  order 
to  determine  the  model's  sensitivity  to  that  parameter. 


Finally,  we  defined  the  local  polarization  contrast: 


Cpot 


Q|  -  Qi 

(0|  +  aJ/2 


(12) 


Because  of  the  orthogonality  of  the  modes  along  the  z  axis, 
we  can  define  the  terms  of  Eq.  (12)  simply:  at  -  2*aml 
and  Ox  “  ZmOatx. 

The  denominator  of  Eq.  (12)  is  proportional  to  the  aver¬ 
age  quantal  absorption  rate  of  two  adjacent  orthogonal 
photoreceptors  (see  Figs.  2  and  3)  under  a  patch  of  inci¬ 
dent  radiation  uniformly  polarized  along  the  x  axis  of  one 
of  them  (and  hence  along  the  y  axis  of  the  other).  The 
numerator  is  proportional  to  the  difference  between  the 
quantal  absorption  rates  of  each  of  these  cells. 


3.  NUMERICAL  RESULTS  OF  THE 
GRADED-INDEX  MODEL 

A.  Comparison  with  the  Step-Index  Model 

The  electric  fields  excited  by  x  and  y  polarized  incident 
fields  have  nearly  the  same  magnitude  (cf.  the  values  of 
the  a„'s  resulting  from  the  integration  of  the  electric  fields 
in  Table  1).  However,  the  magnitudes  are  different,  and 
this  difference  is  quite  possibly  large  enough  to  be  ex¬ 
ploited  by  a  fish’s  nervous  system.  Given  the  parameters 
n(x)|,-3.t»a  “  1.33  and  n,,,  —  1.39,  is  1.6%  when  the 
incident  wavelength  is  620  nm  (Ami  of  the  pigment  is 
621  *  1  nm).7 

In  order  to  determine  what  the  graded-index  profile 
gains  in  polarization  sensitivity  relative  to  a  step-index 
model,  we  performed  a  similar  analysis  for  a  guide  with  a 
step  profile  that  had  the  same  thickness  and  average  re¬ 
fractive  index: 

n  —  1.39  |x|  <  3.5  fim , 

n  -  1.33  |x|  at  3.5  Mm .  (13) 

The  procedure  for  determining  was  essentially  the 
same.  However,  in  this  case  there  is  no  ambiguity  in 
the  decision  regarding  whether  a  mode  is  guided:  the 
0’s  are  found  from  transcendental  equations  that  have  a 
finite  number  of  real  values  at  any  given  wavelength.41 
Just  as  in  the  graded-index-profile  model,  there  are  five 
even-symmetric  modee  propagating  for  each  x-polarized 
and  y -polarized  incident  beam  with  free-space  wavelength 
of  620  nm.  The  value  of  C**  obtained  from  the  analysis  of 
the  step-profile  model  at  this  wavelength  is  0.32%;  the 
graded-index  model  is  thus  five  times  as  sensitive  to  po¬ 
larisation  as  the  step-profile  model. 

B.  Wavelength  Dependence  of  Polsrization  Sensitivity 
The  differential  polarization  sensitivity  of  the  green  sun- 
fish  declines  at  690  nm.*  We  wished  to  explore  the  possi¬ 
bility  that  this  might  be  explained  by  our  model  of  the 
double  cones’  ability  to  guide  light.  Figure  6  is  a  plot  of 
C„i  versus  the  free-space  wavelength  (A)  of  incident  radia¬ 
tion  for  both  the  graded-  and  the  step-index  models.  In 
the  graded-index  model  with  the  parameters  determined 
above,  the  power  available  for  phototransduction  at  rele¬ 
vant  wavelengths  is  always  larger  when  the  incident  light 
is  polarized  along  the y  axis  (C^  it  always  positive);  how¬ 
ever,  there  are  discontinuities  in  C^(A)  that  result  from 
the  loss  of  propagating  modes.  At  discrete  wavelengths, 
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Fig.  &  Polarisation  contrast  as  a  function  of  wavelength  for  one* 
dimarcional  pinbolk-indn-  ind  sttp-indu-profita  modrit  Cfg 
ia  defined  fay  Eq.  (11)  and  is  meant  to  bo  an  indication  of  tbs 
model's  preference  for  one  polarisation  over  the  other.  Positive 
and  negative  values  indicate  a  preference  for  elect  ric/idd*  polar¬ 
ised  along  the  y  axis  and  the  x  axis,  respectively.  For  the  graded- 
indsx  model,  n_  -  1.42  and  x#  -  10.0  pm  in  Eq.  (3).  For  the 
etep- index  model,  n  is  1.39  when  |x|  <  3.6  iu n  and  is  1.33  else. 
Note  that  for  some  wavelengths  is  negative  in  the  case  of  the 
step-index  model,  whereas  iVt  condition  never  obtains  in  the 
graded-index  model.  In  the  graded-indax  model  there  is  a  brief 
interval  near  640  nm  where  y-polarised  radiation  excites  one 
additional  mode  compared  with  x-potarised  radiation.  At  this 
paint  Op*  leaps  to  12.9%;  the  spike  is  clipped  so  that  the  rest  of 
each  function  can  be  shown  at  a  suitable  scale. 


modes  get  cut  off  such  that  fewer  guided  modes  propegate 
at  longer  wavelengths.  However,  the  power  available  for 
phototransduction  from  the  TE  modes  rises  more  quickly 
than  that  for  the  corresponding  TM  inodes.  Just  before 
two  corresponding  TE  and  TM  inodes  are  cut  off,  the  dif¬ 
ference  in  power  that  they  capture  from  the  two  orthogo¬ 
nal  incident  polarizations  becomes  substantial.  After  a 
given  set  of  inodes  is  cut  off,  the  value  of  drops,  since 
the  remaining  aet  of  modes  is  not  so  sensitive  to  polarisa¬ 
tion.  As  the  wavelength  is  increased  further,  C„ i  in¬ 
creases  again  until  the  next-lower  aet  of  modes  becomes 
cut  off.  Another  distinctive  feature  of  the  plot  in  Fig.  6  is 
the  spike  at  the  cutoff  frequency  of  one  of  the  TM  modes. 
The  spike  arises  in  the  graded-index  model  because  each 
TM  mode  ia  cut  off  at  a  slightly  shorter  wavelength  than 
ia  the  corresponding  TE  mode.  Over  a  short  range  of 
wavelengths,  the  number  of  propagating  TE  mode*  is  one 
greater  than  the  number  of  propagating  TM  modes.  For 
this  set  of  parameters  the  spike  is  only  0.4  nm  wide,  how¬ 
ever,  so  we  doubt  that  it  could  significantly  affect  the 
polarization-analyzing  capability  of  the  structure  given 
natural  excitation,  which  ia  typically  broadband.4* 
Contrast  the  waveiangth  dependence  of  the  graded-index 
modal  with  that  of  the  step- index  model  in  the  lower  part 
of  Fig.  6.  The  first  notable  difference  is  that  C^i  is  sig¬ 
nificantly  smaller  for  the  step- index  than  for  the  graded- 
indax  modal.  The  second  feature  of  note  is  that,  in  the 
atap-index  model,  Cp*  reverses  sign  at  some  wavelengths. 
In  parts  of  the  spectrum,  the  energy  extracted  by  the 
ebromophores  from  light  polarised  along  the  x  axis  is 
greater  than  that  extracted  from  light  polarized  along  the 
y  axis.  Notice  also  that  there  art  no  spikaa  in  this  plot, 
bacause  in  the  step-index  slab  the  corresponding  TE  and 
TM  modes  are  cut  off  at  the  same  waveiangth. 


Rowe  rial. 

C  Sensitivity  of  the  Model  to  Variations  in  ia,  and  x. 

We  calculated  the  effect  of  a  higher  average  refractive  in¬ 
dex  by  increasing  the  value  of  the  maximum  refractive 
index  [nM  in  Eq.  (2)  et  seq.].  Given  that  the  index  of  re¬ 
fraction  arises  principally  from  proteins  and  lipids  within 
the  cytoplasm,  is  unlikely  to  be  much  greater  than 
1.6.**  However,  given  thoae  constituents,  it  may  well  be 
higher  than  our  original  assumption  of  1.42.  We  there¬ 
fore  repeated  the  calculations  for  various  values  of  n be¬ 
tween  1.42  end  1.6  end  plotted  them  in  Fig.  7.  For  these 
calculations  we  retained  the  constraint  that  - 

1.33;  so,  as  nw  increased,  x«  decreased  (i  a,  r  fell  off  with 
x  et  e  steeper  rate).  We  continued  to  define  the  cutoff 
point  in  the  same  fashion  (modes  were  cut  off  if  their 
effective  index  wee  less  than  1.33)  and  considered  d  in 
Eqi.  (10)  end  (11)  to  be  7  pm.  Given  these  stipulations,  as 
Ran  increases,  the  number  of  modes  propagating  at  each 
wavelength  increases,  or  at  least  remains  the  same.  The 
relative  differences  between  a. (  end  aml  ere  larger  for 
each  of  the  corresponding  modes,  however,  so  Cpu  ia  larger 
at  each  wavelength  (except  at  the  spikes).  As  more  modes 
ere  propagating,  there  ere  more  inodes  to  get  cut  off,  so 
CpriCA)  has  more  jumps  in  it  vdien  R.U  is  larger.  However, 
since  individual  modes  do  not  contribute  so  much  to  the 
total  absorbed  energy,  the  relative  magnitudes  of  the 
jumps  are  not  as  large  for  the  larger  values  of  Ran- 
Because  the  value  of  the  thickness  that  we  used  [d  in 
Eqe.  (10)  and  (11)]  was  baaed  on  an  upper  limit  for  the 
cross  section  of  an  inner  segment,  we  wanted  to  deter¬ 
mine  how  polarization  contrast  would  behave  in  e  smeller 
structure.  Another  reason  for  determining  this  was  that 
we  believed  that  it  might  be  more  appropriate  to  think  of 
the  ellipsoid  rather  than  the  entire  inner  segment  as  the 
guiding  structure.  Since  the  boundaries  of  the  ellipsoid 
are  contained  within  the  boundaries  of  the  cytoplasmic 
envelope,  the  ellipsoid  always  presents  e  smaller  cross  sec¬ 
tion  than  does  the  inner  segment  in  a  given  section.  We 
explored  the  effect  of  a  smaller  size  by  fixing  at  1.42 
but  choosing  a  smaller  value  of  x  at  which  to  pin  the  index 
at  1.33.  The  narrowest  portion  of  the  inner  segment  is 
the  scleral  end,  just  before  the  junction  with  the  outer  seg- 


Fig.  7.  Polarization  contrast  as  a  function  of  wavelength  for  dif¬ 
ferent  value*  of  maximum  index  in  the  parabolic-profile  model 
In  each  case,  the  index  profile  satisfies  Eq.  (3)  with  x«  set  such 
that  n  -  1.33.  The  values  of  n^.  are  depicted  at  the  right  of  the 
curve*.  Aa  in  Fig.  6,  the  number  at  the  top  of  each  spike  repre¬ 
sent*  that  value  to  which  CN  jump*  at  that  wavelength  for  the 
cropped  curves. 
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Fig.  8.  Polarisation  contrast  as  a  function  of  wavelength  for  dif- 
foront  values  of  thickness  in  the  parabolic -profile  modal.  For 
ease  of  viewing,  the  curves  are  plotted  on  two  different  sots  of 
axes  (A  and  B).  For  each  curve,  the  index  profile  satisfies  Eq.  (3) 
with  iimm  m  1-42.  The  vahiss  of  |*|  where  n  -  1.33  art  consid¬ 
ered  the  walls  of  the  guide  such  that  the  thickness  is  twice  thst 
value.  Depicted  at  the  right  of  each  curve  is  this  thick¬ 
ness.  As  in  Figs.  6  and  7,  the  number  at  the  top  of  each  spike 
represents  the  value  to  which  Jumps  at  that  wavelength  for 
the  cropped  curves.  For  the  curve  with  thickness  equal  to  2  ju a, 
the  spike  is  4  am  wide  but  is  roughly  flat,  with  C,*  -  384%.  All 
other  spikes  are  less  than  1  nm  wide. 


meat  (pm  Fig.  1).  At  this  region  of  the  double  cone,  the 
fOtpuid  may  be  as  small  2  or  3  pm  wide  along  the  short¬ 
est  axis.  We  therefore  chose  to  perform  the  calculations 
again  tar  vahiss  of  d  between  2  and  7  pm.  The  results  of 


these  calculations  are  displayed  in  Fig.  8.  As  expected, 
fewer  modes  propegate  es  the  thickness  of  the  wsveguide 
approaches  the  size  of  the  wavelength  of  incident  light. 
When  fewer  mode*  are  propagating,  the  loas  of  e  single 
mode  has  a  much  greater  effect  on  thus  general  state¬ 
ments  about  the  usefulness  of  smaller  structures  are  some¬ 
what  more  problematic  than  statement#  about  changes  in 
n—..  Decreasing  d  also  makes  the  spikes  tend  to  become 
higher  and  wider.  The  spike  in  the  case  where  d  -  2  pm 
jumps  to  nearly  40%  end  is  4  nm  wide.  Since  this  spike 
is  near  the  wavelength  of  maximum  absorption  for  the 
double-cone  pigment,  it  might  have  a  large  effect  on 
the  polarization-analysing  ability  of  the  photoreceptor. 
Therefore  in  the  next  paragraph  we  develop  another  quan¬ 
tity  that  incorporates  a  consideration  of  the  double  cones’ 
sensitivity  across  wavelength. 

Broadband  polarization  contrast.  To  estimate  the  ef¬ 
fect  of  broadband  polarized  light  on  the  relative  polariza¬ 
tion  difference  in  outputs  between  orthogonal  double  cones 
in  the  sunfish  retina,  we  define  the  quantity 


■ 
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where  A  is  the  free-space  wavelength  in  nanometers  end 
SKA)  it  the  action  spectrum  of  the  double-cone  pigment 
(data  collected  for  Ref.  7,  and  courtesy  of  Allen  Dearry).45 
The  Pim  used  in  the  determination  of  the  a’s  was  assumed 
to  be  a  constant  across  wavelength.  The  values  of 
thus  calculated  can  be  found  in  Tables  2  and  3. 


Table  S.  Integration  of  Differential  Polarization 
Sensitivity:  d  Fixed  at  7  pm’ 


Ami 

^,<%)* 

1.42 

1.480 

1.46 

1.915 

1.48 

2.341 

1.81 

2.737 

1.64 

3.148 

1.57 

3.557 

1.60 

3.966 

*As  in  tho  plots  of  Fig.  7. 

TLsi  is  dsfinsd  in  Kq.  (14).  In  mm 

nee  it  is  a  measure  of  tho  contrast 

gonorstod  by  s  psir  of  orthogonal  double  cooas  stimulated  by  a  broadband 

beam  of  light,  linearly  polarised  parallel  to  tho  long  axis  of  tho  cross  section 
of  one  of  those  cones. 

Table  3.  Integration  of  Differential  Polarisation 


Sensitivity: 

Ran  Fixed  at  1.42* 

d  {pm) 

C*<%)* 

2 

2.742 

3 

1.386 

4 

1.647 

6 

1.527 

6 

1.454 

7 

1.460 

*  A*  in  tbs  plots  of  Fig.  8. 
"tft  is  dofinsd  in  Eq.  (14). 
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4.  DISCUSSION  AND  CONCLUSIONS 

A.  Resolution  of  Throe  Puzzles 

Wb  found  that  a  graded-index  model  of  the  sunfish  double¬ 
cone  inner  segment  predicts  the  inner  segment  to  guide 
light  in  a  polarisation-dependent  manner.  We  hypothe¬ 
sise  that  such  polarisation-dependent  guiding,  or  geomet¬ 
ric  birefringence,  could  provide  the  physical  basis  for 
differential  polarization  sensitivity  in  many  vertebrates. 
This  hypothesis  links  and  could  resolve  three  long-standing 
pussies:  (1)  the  puzzle  of  double  cones,  (2)  the  puzzle  of 
double-cone  mosaics,  and  (3)  the  puzzle  of  vertebrate  po¬ 
larisation  sensitivity. 

The  puzzle  of  double  conet.  In  discussing  double  cones, 
Walls  states  (p.  58  of  Ref.  9)  that  “the  puzzle  [that  double 
cones]  present  is  particularly  irritating  because  they  are 
to  very  widespread  among  vertebrates,”  and  that  “They 
have  unfortunately  not  greatly  interested  visual  physiolo¬ 
gists  since  the  latter  have  their  attention  focussed  upon 
the  human  retina,  in  which  double  cones  are  lacking.”  In 
the  60  years  intervening  since  Walls  wrote  these  state¬ 
ments,  there  has  been  no  hypothesis  that  attributes  to 
double  cones  a  utility  distinct  from  that  of  single  photo¬ 
receptors.  The  hypothesis  that  double  cones  function  to 
take  advantage  of  polarization,  a  fundamental  property  of 
light,  is  appealing.  If  this  hypothesis  is  correct,  it  sug¬ 
gests  that  many  retinas  take  advantage  of  a  wealth  of 
information  contained  in  the  polarization  of  light  consti¬ 
tuting  a  visual  scene,  information  of  which  we  are  only 
beginning  to  become  aware.44*41 

The  puttie  of  double-cone  moeaict.  The  regularity  of 
the  double-cone  mosaics  in  many  vertebrates  is  striking 
(see  Refs.  15-18  and  also  Fig.  3).  It  has  been  hypothe¬ 
sised  that  these  mosaics  form  the  basis  of  a  movement- 
detection  system.4*  However,  motion  detection  by  retinas 
that  have  no  double  cones  is  well  established.4*  More  im¬ 
portant,  given  our  conclusion  that  double  cones  guide  light 
as  unitary  structures,  the  double-cone  mosaic  would  be 
lese  suitable  for  motion  detection  than  a  similar  retinal 
mosaic  with  each  double  cone  replaced  by  two  single  cones. 
Motion  detection  requires  the  comparison  of  signals  from 
neighboring  receptors.  Analysis  of  Figs.  1  and  2  leads  to 
the  conclusion  that  light  impinging  upon  one  of  the  mem¬ 
bers  of  a  double-cone  pair  will  be  guided  into  the  outer  seg¬ 
ments  of  both  members,  so  a  comparison  of  the  responses 
of  the  members  of  a  double  cone  would  be  less  informative 
about  motion  than  comparison  of  the  responses  of  single 
cones  with  the  same  packing  density.  On  the  contrary, 
the  puzzle  of  the  orderly  mosaics  of  double  cones  may  come 
into  sharper  focus  in  light  of  the  hypothesis  that  double 
cones  are  geometric  birefringent  polarization  analyzers. 
Under  this  hypothesis  the  double-cone  mosaics  could  serve 
as  the  input  for  the  calculation  of  a  polarisation-contrast 
or  -difference  image,  analogous  to  a  color-contrast  image 
computed  by  color-opponent  cells  (see  below). 

Although  the  regularity  of  teleost  double- cone  mosaics 
has  suggested  to  others  a  function  in  polarisation,1*  the 
general  consequences  of  this  idea  seem  not  to  have  been 
pursued.  One  such  consequence  is  that,  in  addition  to  elu¬ 
cidating  a  role  for  the  square  mosaics  such  as  in  sunfish, 
the  hypothesis  also  provides  an  explanation  for  the  un¬ 
usual  mosaics  described  by  Lyall  in  the  pike*1  and  by  Walls 
in  the  scorpion  fish.*  In  these  animals  the  long  axes  of 


adjacent  double  cones  are  offset  from  each  other  by  60*. 
Under  the  hypothesis  that  double  cones  are  birefringent 
polarization  analyzers,  the  mosaic  of  the  pike’s  retina  con¬ 
tains  a  three-channel  polarization-analyzing  system  analo¬ 
gous  to  a  trichromatic  color-vision  system  (see  Ref.  19). 

The  puzzle  of  differential  polarization  teniitivity. 
Waterman’s  discovery  nearly  35  years  ago*9  and  subse¬ 
quent  research  demonstrating  vertebrate  polarization 
sensitivity1"4  have  given  us  the  puzzle  of  a  capacity,  polar¬ 
ization  sensitivity,  with  no  established  biophysical  mecha¬ 
nism  (there  is,  however,  one  notable  exception  in  which  the 
mechanism  has  been  established).69  It  has  been  shown 
through  a  variety  of  techniques  that  the  photoreceptors  of 
primates  dynamically  orient  their  longitudinal  axes  to¬ 
ward  the  exit  pupil  of  the  eye.6441  This  orientation  serves 
to  maximize  the  light  guided  to  the  outer  segments,  since 
off-axis  light  is  ineffectively  trapped  by  photoreceptor 
inner  segments,  a  phenomenon  that  explains  the  Stiles- 
Crawford  effect  of  the  first  kind.**  However,  such  orien¬ 
tation  of  the  photoreceptors  toward  the  exit  pupil  also 
minimizes  the  photoreceptors’  potential  use  of  the  dichro- 
ism  inherent  in  the  photopigments  and  outer  segment 
disks  for  polarization  discrimination.*0  This  is  so  became 
the  E-vector  of  axially  guided  light  is  parallel  primarily  to 
the  planes  of  the  outer  segments  disks,  in  the  same  planes 
as  the  axes  of  the  chromophores.  Although  it  has  yet  to 
be  established  that  the  longitudinal  axes  of  fish  photo¬ 
receptors  are  oriented  toward  their  entrance  pupils,  it 
aeems  likely,  on  the  grounds  of  optimization  of  visual  sen¬ 
sitivity,  that  they  are.  Certainly,  the  establishment  of 
waveguiding  behavior  in  teleost  cones19  and  the  measured 
directional  sensitivity  of  turtle  cones*1  suggest  that  photo¬ 
receptor  orientation  in  these  animals  will  follow  a  pattern 
similar  to  that  established  in  primates. 

Nonetheless,  some  researchers  have  hypothesized  that 
the  dichroism  of  the  photopigments  in  zitu  could  serve  for 
polarization  detection  by  the  capture  of  light  that  is  not 
propagating  radially  through  photoreceptors.  For  ex¬ 
ample,  it  has  been  proposed  that  polarization  detection 
might  occur  in  the  embryonic  fissure,  a  prominent  struc¬ 
ture  in  many  fish.*1  In  this  region  of  the  retina,  light 
may  strike  photoreceptors  from  the  sides,  allowing  the 
photoreceptors  to  make  use  of  the  dichroism  conferred  by 
the  stacking  of  the  disks  and  the  transition  moments  of 
the  photopigments.  However,  differential  polarization 
sensitivity  was  demonstrated  in  the  sunfish  when  the 
stimulus  was  restricted  to  a  region  of  the  retina  remote 
from  this  fissure.*  In  the  remainder  of  this  discussion  we 
explore  ramifications  of  the  hypotheses  proposed  here, 
that  the  inner  segments  of  the  photoreceptors  in  the  green 
sunfish  have  a  graded  refractive-index  profile  and  that  the 
three-dimensional  shape  of  this  profile  gives  the  photo¬ 
receptors  sufficient  geometric  birefringence  for  that  bi¬ 
refringence  to  serve  as  a  biophysical  basis  for  polarization 
sensitivity. 

B.  Graded-Index-Profile  Waveguide  Model  of  the 
Double-Cone  Inner  Segment 

Exelution  of  the  Hep-index  model.  With  the  size  and  re¬ 
fractive  indices  expected  for  the  green  sunfish  double-cone 
inner  segment,  the  polarisation  contrast  of  the  step-index 
model  was  approximately  0.3%  for  wavelengths  ranging 
from  570  to  620  nm.  It  is  possible  that  an  animal  could 
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take  idnotifi  of  such  small  differences  between  the 
quantal  catch  rates  of  its  different  cone  classes:  in  rhesus 
macaques  the  increment  threshold  for  a  email  monochro¬ 
matic  flash  on  a  30,000-Td  white  adapting  background 
is  reached  when  the  chromatic  contrast  between  the 
long*  versus  middle-wavelength-sensitive  cones  is  0.28%  or 
0.6*.“  For  humans  who  are  light  adapted  to  a  less- 
intense  field,  however,  the  just-detectable  chromatic  con¬ 
trast  between  the  long-  and  middle- wavelength-sensitive 
cones  is  approximately  4%.“  Thus  we  doubt  that  a  sun- 
fish  could  And  a  maximal  cone-contrast  difference  of 
0.3%  in  its  double  cones  generally  useful,  given  that  re¬ 
flected  light  from  natural  targets  is  unlikely  to  be  per¬ 
fectly  polarised. 

An  additional  argument  against  the  step-index  guide 
can  be  baaed  on  the  sign  reversals  of  Cp*  (see  Fig.  6). 
Such  sign  reversals  would  decrease  the  utility  of  a  double¬ 
cone  differencing  mechanism  in  real-world  conditions,  in 
which  radiation  tends  to  be  broadband  and  polarization  is 
a  weak  function  of  wavelength.  Within  a  photoreceptor, 
polarization  sensitivity  at  specific  wavelengths  will  be 
offset  fay  any  opposite-signed  polarization  sensitivity  at 
other  wavelengths. 

Basts  for  geometric  birefringence.  The  evidence  for  a 
graded-index  profile  in  sunfish  double  cones  suggested 
to  us  that,  despite  their  relatively  large  size,  the  double 
cones  might  be  polarization  sensitive  as  a  result  of  geomet¬ 
ric  birefringence.  When  a  dielectric  waveguide  is  com¬ 
posed  of  a  nonchiral  and  isotropic  material,  its  polarization 
properties  depend  only  on  its  shape.  However,  the  shape 
of  a  waveguide  is  generally  important  for  polarization- 
dependent  behavior  only  when  the  dimensions  of  the  guide 
are  comparable  with  the  wavelength  of  incident  radiation. 
The  narrowest  dimensions  of  the  sunfish  double-cone 
cross  sections  are  typically  5-10  times  the  wavelength  of 
physiologically  relevant  light.  If  the  index  is  graded  in 
the  manner  that  we  propose,  however,  then  the  guided 
power  is  concentrated  closer  to  the  central  axis  of  the 
double  cone  than  it  would  be  if  the  index  were  homoge¬ 
neous.  Thus  a  graded-index  waveguide  is  effectively 
smaller  than  a  step-profile  guide  with  the  same  outer 
dimensions  and  average  refractive  index,  and  therefore  its 
polarization  contrast  is  greater. 

Abaence  of  extracellular  cladding.  A  striking  feature 
apparent  in  the  electron  micrographs  of  the  photoreceptors 
is  that,  near  the  base  of  the  inner  segment  (the  vitreal  end), 
there  is  an  absence  of  any  significant  extracellular  space. 
In  other  words,  it  appears  that  the  inner  segments  occupy 
virtually  the  entire  retinal  cross  section  (see  Fig.  1  and 
Figs.  2b-2e).  In  most  modern  analyses  of  photoreceptors 
as  waveguides,  the  index  profile  is  assumed  to  be  a  step 
that  implicitly  requires  an  extracellular  cladding  (of  thick¬ 
ness  comparable  with  at  least  the  relevant  wavelengths, 
a  condition  that  clearly  is  not  met;  see  the  620- nm  scale 
bar  in  Fig.  3a)  in  order  to  sustain  bound  modes  easily.  In 
our  mathematical  analysis  of  a  waveguide  model  with  a 
graded-index  profile,  we  have  also  considered  the  extracel¬ 
lular  region  as  cladding  with  a  refractive  index  of  1.33. 
However,  as  can  be  seen  from  Fig.  6,  the  excited  bound 
modes  are  not  very  sensitive  to  the  presence  of  this  extra¬ 
cellular  cladding,  since  they  die  out  rapidly  in  the  region 
beyond  die  waveguide’s  walls.  This  is  due  mainly  to  the 
fact  that  in  this  graded-index -profile  model  the  power  of 


a  mode  is  concentrated  near  the  guide  axis  where  the 
index  of  refraction  is  the  highest.  This  might  suggest  a 
new  approach  to  the  modeling  of  photoreceptor  wave¬ 
guides  and  has  ramifications  for  the  analysis  of  visual 
acuity,  our  next  topic. 

Possible  role  in  vieual  acuity.  A  factor  that  potentially 
limits  visual  acuity  in  any  species  is  the  interphoto¬ 
receptor  spacing:  for  maximal  acuity,  photoreceptors  of 
any  fixed  size  should  be  tightly  packed.  Indeed,  the  spac¬ 
ing  between  photoreceptors  is  thought  to  be  the  major  fac¬ 
tor  limiting  acuity  in  fish  visual  systems **-“  and  perhaps 
in  a  variety  of  other  animals,  including  humans  (under 
some  conditions).**’0  Tight  packing,  however,  necessarily 
increases  energy  coupling  between  receptors,  reducing 
contrast  between  adjacent  cells  and  hence  reducing  visual 
acuity.  Perhaps  the  greatest  advantage  of  a  graded-index 
profile  is  the  concentration  of  power  near  the  central  guide 
axis.  This  power  concentration  may  serve  to  isolate  photo¬ 
receptors  optically  from  one  another,  while  permitting 
maximal  packing.  A  graded  refractive-index-profile 
guide  draws  light  toward  its  center,  reducing  coupling  and 
permitting  higher  acuity  for  a  given  receptor  size  and 
spacing.  This  advantage  of  a  graded-index  profile  does 
not  depend  on  the  noncircular  cross  sections  of  the  double 
cones.  It  would  also  be  conferred  on  the  single  cones, 
which,  as  seen  in  Fig.  2,  also  exhibit  a  differential  electron 
density  across  their  ellipsoids. 

Reduction  of  backreflectiona  by  inner  tegmenta.  In 
both  single  and  double  cones,  the  electron  density  appar¬ 
ently  rises  and  then  levels  off  in  the  progression  from  the 
base  of  the  ellipsoid  to  the  base  of  the  outer  segments 
(see  Fig.  1).  An  advantage  conferred  by  such  z-axis  gra¬ 
dation  is  a  reduction  in  the  reflection  of  axially  propagat¬ 
ing  light.  Sharp  changes  in  refractive  index  represent 
impedance  boundaries  that  produce  backreflectiona.  A 
gradual  increase  in  refractive  index  reduces  these  back- 
reflections.  Many  insect  ommatidia  are  adapted  to  take 
advantage  of  this  effect;  dimples  on  the  surfaces  of  their 
corneas  serve  as  antireflection  coatings  by  reducing  the 
rate  at  which  the  refractive  index  changes  as  light  travels 
into  the  animals’  eyes.’1,7*  Sunfish  photoreceptors  ap¬ 
pear  to  take  advantage  of  the  same  principle  in  their  photo¬ 
receptor  inner  segments. 

Combined  with  the  gradation  in  index  along  the  z  axis, 
the  x  and  y  axes’  gradation  may  reduce  backscattering  of 
nonparaxial  light.  Fish  have  large  pupils  and  nearly 
spherical  lenses.  The  numerical  aperture  of  the  goldfish 
eye  is  approximately  0.44,“  allowing  the  photoreceptors  to 
collect  light  incident  up  to  ±19*  off  of  their  longitudinal 
axes.71  For  reduction  of  the  reflection  at  an  interface,  it 
is  desirable  to  grade  the  index  slowly  in  the  direction  along 
which  light  is  propagating.  Thus  the  three-dimensional 
index  profile  of  sunfish  cones  may  minimize  the  reflec¬ 
tion  of  light  over  a  large  entrance  aperture  while  also 
concentrating  the  guided  light  near  the  center  of  the  guid¬ 
ing  structure. 

C.  Role  of  the  Double-Cone  Mosaic:  A  Basis  for 
Polarization  Opponency 

Figure  9  diagrams  a  neural  mechanism  by  which  large 
differential  polarization  sensitivity  can  arise  from  two 
populations  of  cells  whose  sensitivity  is  a  weakly  modu¬ 
lated  function  of  polarization.  The  two  sinusoids  (dashed 
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Fig.  ft  Opponency  between  the  output*  of  two  orientation  cla*M* 
of  double  cute  can  drastically  affect  the  differential  polarization 
sensitivity  of  the  whole  animal.  The  output  of  a  double  cone 
might  be  only  weakly  modulated  by  polarization,  a*  indicated  by 
the  two  dashed  curve*  at  the  top.  The  shading  in  the  schematic 
double  cones  indicates  which  orientation  class  should  be  maxi¬ 
mally  stimulated  at  those  angles.  The  large  solid  curve  in  the 
figure  represents  the  unsigned  difference  between  the  outputs  of 
the  two  sets  of  double  cones.  Each  of  the  three  curves  r  normal¬ 
ized  according  to  its  maximum  output. 

curves)  indicate  the  response  of  individual  double  cones  as 
a  function  of  polarization  angle  to  linearly  polarized  light 
of  constant  intensity,  as  computed  for  Table  1.  The  sche¬ 
matic  images  of  a  double-cone  tetrad  indicate  which  cells 
give  the  largest  response  at  their  respective  maxima. 
The  solid  curve  represente  the  absolute  value  of  the  differ¬ 
ence  between  the  outputs  of  the  receptors,  normalized  to 
the  maximum  difference.  An  inner  retinal  neuron  that 
receives  output  from  the  receptors  in  an  opponent  fashion 
would  exhibit  responses  with  this  magnitude  profile,  and 
measurements  of  an  animal’s  sensitivity  could  be  ex¬ 
pected  to  follow  such  a  pattern.  Evidence  for  polarisation 
opponency  in  the  sunfish  retina  has  been  obtained  in 
behavioral  experiments.*  Furthermore,  although  pure 
polarization-opponent  cells74  have  yet  to  be  found  in  elec- 
trophysiological  experiments  in  a  vertebrate,  they  have 
been  predicted  to  exist  in  arthropods  and  have  been  found 
in  the  cricket.71  We  predict  that  similarly  responsive 
cells  will  be  found  in  the  retina  of  a  vertebrate. 
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